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1. Executive Summary 

Herein we summarize the major results of the five-year MURI, "Prevention of 
Corrosion in Structural Aluminum Alloys." This program focused on examining 
mechanisms of copper redistribution during corrosion and protection of structural 
aluminum alloys such as 2024. Our major accomplishments are summarized in the 
bullets below. 

• The development of a new electrochemical technique based on Pb underpotential 
deposition for the quantitative determination of copper redistribution during corrosion 
and coating pretreatments of copper containing aluminum alloys. 

• The development of synthetic (model) 2024-T3 A1 alloys that can be used to simulate 
the corrosion behavior of the real 2024-T3 alloy. Synthetic 2024-T3 is made using 
conventional lithography. 

• The determination that our “synthetic A1 2024-T3” behaves as efficiently as a planar 
copper electrode for oxygen. Therefore, the synthetic alloy serves as an excellent 
model system for the cathodic reactions during corrosion of real 2024-TS. 

• Characterization of the dealloying and corrosion behavior of the AECuMg 
intermetallic S-phase which is the dominant second phase particle present in A1 Alloy 
2024-T3. 

• The development of a new technique which uses the tip of a scanning tunneling 
microscope (STM) as a local pH probe at nanometer length scales. The purpose of 
this study was to determine the relative activity of the various particles on the surface 
of A1 alloy 2024 for oxygen reduction. Our results indicated that the S-phase is 




2 


dominant with respect to controlling oxygen reduction on the surface of corroding A1 
2024-T3. 

• The development of a kinetic Monte Carlo code that simulates dealloying current- 
voltage behavior and morphology in any single-phase alloy system. This code was 
used to examine dealloying of copper in aluminum solid solution and dealloying of S- 
phase. 

• Identification that during corrosion of A1 Alloy 2024-T3, the S-phase is first anodic 
(i.e., it undergoes dealloying) and then turns in to an active cathode for oxygen 
reduction that supports corrosion. 

• Determination of the diffusion - limited cathodic current density from oxygen 
reduction on the surface of A1 alloy 2024 T-3 (20 (lAcm" 2 ). 

• Determination that our “synthetic A1 2024-T3” behaves as efficiently as a planar 
copper electrode for oxygen. Therefore, the synthetic alloy serves as an excellent 
model system for the cathodic reactions during corrosion of real 2024-TS. 

• Results of Rotation Ring-Disk electrode experiments demonstrated the operation of a 
dissolution/back-plating copper redistribution mechanism during corrosion of A1 
2024-T3. This process is triggered by very small copper particles of order tens of 
nanometers in size that become mechanically disconnected from de-alloyed S-phase 
intermetallic particles. As long as the remnant de-alloyed copper sponge stays 
mechanically and electrically connected to the substrate no copper dissolution can 
take place. In this case, we can disregard particle size effects on dissolution since the 
corrosion potential is almost 1 V less than the reversible potential of copper. 

• Results of RRDE experiments demonstrated that the most significant source of 
redistributed copper during corrosion and pretreatment of A1 alloy 2024 comes from 
matrix dealloying. This is in agreement with results of atomic - scale Kinetic Monte 
Carlo simulations. The other source of redistributed copper is mechanical 
disconnection of copper (from dealloyed S-phase) that dissolves under open circuit in 
the electrolyte and subsequently electrochemically deposits on the surface of A1 alloy 
2024. 

• Realization that the low limiting net cathodic current density (~ 0.5 |iAcm' 2 ) 
delivered by Al-Cu solid solutions is due to the superimposition of alkaline anodic 
dissolution on the cathodic oxygen reduction reaction, a subtle example of so-called 
“cathodic corrosion”. The system is self-regulating in that the net current is just 
enough to poise the surface pH at the value where the open-circuit potential of the 
alloy is equal to the applied or naturally prevailing potential. 

• Finally, the differences in behavior for the surface coverage of Cu between 2024-T3 
and the synthetic 2024 alloy is attributed to the operation of a dissolution/back-plating 
mechanism triggered by mechanical disconnection of Cu particles from de-alloyed 
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porous S-phase. Comparisons of the behavior of these two alloys allow us to 
conclude that the matrix de-alloying and dissolution/back-plating contribute about 
equally to the copper redistribution process during corrosion of Aluminum Alloy 
2024-T3 in 0.5 M NaCl. 
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5. Major Accomplishments 

• Nature publication related to fundamental aspects of de-alloying. 

• The result of this MURI investigation has been the identification of a new corrosion 
mechanism, i.e., matrix de-alloying. This type of process allows for renewal of 
cathodic sites during corrosion of not just Al alloys, but for all materials. We believe 
that this is the significant accomplishment of this 5-year study. 

6. Honors, Awards, Patents 

There have been no honors, awards or patents to date resulting from this work. 
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Cu NANOPARTICLE FORMATION: COPPER REDISTRIBUTION DURING 
NaCl SOLUTION CORROSION OF Al-Cu-Mg ALLOYS 

R. G. Ford, R. W. Carpenter, & K. SieradzkF 

Science & Engineering of Materials Program Center for Solid State Sc ™ ce ’ 8 * 7 '^ eChaiUCal & 
Aerospace Engineering Department, Arizona State University, Tempe, AZ 85287-1704 
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Fig 1. A reaction zone formed on homogenized and aged (20 hrs @ 190°C) for S’ in 2024 alloy, (a) 
Reaction zone and edge of foil to the right. Fe-containing dispersoid marked A and Mil-containing 
dispersoid marked B, both inside the reaction zone, (b) Higher magnification image of the boxed area 
showing smaller nanospheres following <112> matrix orientation of original S* morphology, at 100 
keV. 

Fig 2. Reaction zone at edge of S crystal foil after 5 seconds in 0.5 M NaCl aqueous solution, at 100 
keV. 

Fig 3. Cu-rich nanosphere distribution on the surface of a homogenized and quenched (24 hrs (s' 
495°C) 2024 alloy specimen after 30 sec in 0.5 M NaCl Uq) solution, imaged at 200 keV. 





Microsc. Microanaj. 4 (Suppl 2: Proceedings) 
© Microscopy Society of America 1998 
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Copper Redistribution during Corrosion of Aluminum Alloys 

N. Dimitrov, J. A. Mann, and K. Sieradzki* 

Department of Mechanical and Aerospace Engineering and Science of Engineering Materials Program, Arizona State University, 
Tempe, Arizona 85287-6106, USA 


We describe a new quantitative electrochemical technique for measuring the copper redistribution that occurs during corrosion and 
pretreatment of high-strength copper-containing aluminum alloys. Underpotential deposition (upd) of Pb on Cu was used to assay 
the surface coverage of Cu. Model experiments using samples of known Al-Cu ratios were used to evaluate the accuracy of the 
technique and develop a calibration curve. The results demonstrated the accuracy of the method to within 2%. The elemental cop¬ 
per remaining on the surface of A1 2024-T3 samples after various pretreatment/corrosion procedures was determined using the upd 
technique. A simple treatment involving open-circuit etching of the alloy surface in 3 M HN0 3 was found to be most efficient for 
the removal of Cu from the A1 alloy surface. In accord with a recent observation of Bucheit et al., 1 convection was found to enhance 
copper redistribution on an A1 2024-T3 sample in 0.5 M NaCl solution. 

© 1999 The Electrochemical Society. S0013-4651(98)04-027-0. All rights reserved. 
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High strength is achieved in aluminum aircraft alloys by the pres¬ 
ence of Cu as an alloying element. However, copper is particularly 
deleterious to the corrosion behavior of A1 alloys owing to galvanic 
couples set up between copper-rich intermetallic inclusions and the 
adjacent copper-depleted matrix. This makes copper-containing AI 
alloys the most difficult to protect with conversion coatings. Owing 
primarily to these arguments, it has proven particularly difficult to 
identify a suitable environmentally acceptable replacement for 
chromium-containing conversion coatings for the aluminum alloy 
2024-T3. The corrosion of this alloy is thought to be enhanced by the 
cathodic activity of intermetallic particles. These particles serve as 
preferential sites for the reduction of oxygen, since the conductivity 
of the mixed oxide film covering them is higher than that of the pas¬ 
sive film on elemental aluminum. A key issue in the coating and cor¬ 
rosion protection of Al-Cu-Mg alloys is the redistribution of copper 
that occurs during pretreatment and subsequent corrosion. A signif¬ 
icant accumulation of copper on the alloy surface was already shown 
after mechanical polishing of Al 2024-T3 in aqueous diamond sus¬ 
pension 6 and after chemical polishing in a mixture of concentrated 
sulfuric, phosphoric, and nitric acids. 

Enrichment of copper on the alloy surface is believed to result 
from dealloying of intermetallic particles, in particular S-phase parti¬ 
cles (Al 2 CuMg) which account for —60% of the particle population 
on the alloy surface. The dealloying process results in nanoporous 
Cu-rich remnants that induce pitting around their periphery. 

The most common approach for examining copper redistribution 
on the surface of a high strength Al-Cu-Mg alloys would be to use 
ultrahigh vacuum surface techniques such as Rutherford backscat- 
tering 7,8 or Auger electron spectroscopy (AES ). 9,10 However, the 
copper which redistributes assumes several morphologies which are 
not particularly suited to conventional surface analysis approaches. 
Dealloying produces nanoporous Cu structures and most of this area 
is inaccessible to standard surface analysis. AES would sense the 
projected area of the porous Cu but would give no hint regarding the 
true physical area comprising the pore structure. It is the physical 
area taken up by the elemental Cu that we would like to measure, 
since this determines the cathode area available for oxygen reduc¬ 
tion. Additionally, redistributed monolayer (ML) levels of Cu will 
likely oxidize during transport to a vacuum system unless heroic 
measures are employed. These arguments suggest that an electro¬ 
chemical technique, sensitive to the same physical area of elemental 
Cu that dissolved oxygen would see, would be the preferred method 
of assaying the Cu surface coverage. Such a technique could be 
based upon an electrochemical reaction which is sensitive to copper 
but not to the aluminum or other constituents of the alloy. 


* Electrochemical Society Active Member. 


An appropriate reaction for this purpose is Pb underpotential 
deposition on copper resulting in the formation of one monolayer of 
Pb adatoms on the Cu surface. Pb upd on Cu was investigated on Cu 
polycrystalline electrodes 11 and on monocrystalline Cu (11 1) 12-16 and 
Cu (001 ) 17 faces. It has been shown that Pb forms one epitaxial mono- 
layer on the copper surface 16,17 regardless of the substrate crystallo¬ 
graphic orientation. Twin-electrode thin layer voltammetry suggests 
an ideal charge stoichiometry, i.e., Z = 2 for the system Pb 2+ /Cu 
(111), CIO 4 . 13 Charge densities in the range of 300-335 p-C cm -2 
were reported for the Pb upd ML formed on the Cu (111) elec¬ 
trode 13 ’^ 6,18 and a charge density of 260 (jlC cm -2 for the one formed 
on the Cu (001) electrode . 17 The kinetics of the Pb underpotential 
deposition on Cu (111) were found to depend strongly on the presence 
of a removable oxygen species on the substrate surface . 13,15 A pro¬ 
nounced irreversibility of the process in a solution containing CIO 4 
was reported . 16 The difference in the potentials of the adsorption and 
stripping peak, —150 mV, was attributed to adsorption/desorption of 
oxygen-containing species . 16 Even traces of Cl", CH 3 COO", and 
other specifically adsorbed ions were found to reduce the peak poten¬ 
tial difference and to enhance the voltammetric peak shape . 12,15,16 
The significant influence of the specifically adsorbed ions on the upd 
process kinetics is related mainly to a competition between those 
anions and the oxygen species originally existing on the copper sur¬ 
face. A nucleation-like behavior in terms of a “negative 2D nucleation 
and growth process” within a preformed surface layer consisting of 
specifically adsorbed ions or oxygen-containing species was suggest¬ 
ed for the electrochemical system considered . 15 

Experimental 

We used a model Al-Cu system to access the quality of a quanti¬ 
tative assay of Cu based on the charge associated with the Pb upd 
process. Strips of commercially pure polycrystalline Al were used 
along with strips of commercially pure polycrystalline Cu as a work¬ 
ing electrode in the experiment. Two samples were prepared, each 
with 1 cm 2 Al and 1 cm 2 Cu exposed. The area exposed was polished 
with Buehler Carbimet paper to 15 jxm. On the first sample the 
exposed Al was gradually masked, 0.25 cm 2 at a time, with nail pol¬ 
ish so that results could be obtained for percentages of Cu between 
50 and 100%. On the other sample the exposed Cu was masked in 
the same way, giving results between 50 and 0% Cu. The nail polish 
was stable and showed no signs of decomposition in the electrolytes 
used over the time spans relevant to these experiments. The samples 
were treated for 15 min in H 2 0 2 and dipped for 10 s in 1 M HN0 3 . 
A lead wire etched prior to the experiment in dilute nitric acid was 
used as a reference electrode. A platinum wire coil (area of 3 cm 2 ), 
hydrogen flame annealed prior the experiment, served as a counter 
electrode. The solution used for this electrochemical work consisted 
of KT 1 M NaCI0 4 + 10 -3 M Pb(C10 4 ) 2 + 1(T 4 M HC10 4 + 1(T 3 
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M CH 3 COONa. All the solutions were made with high purity grade 
chemicals and >18 Mfl Nanopure water. A BAS CV-27 model 
potentiostat was used for the experiment and a Nicolet 310 digital 
oscilloscope was used to store data. In each experiment, the poten¬ 
tial was swept from —20 to 300 mV (Pb/Pb 2+ ) at a sweep rate of 
20mV s _1 . All potentials are reported vs. the Pb/Pb 2+ pseudo-refer¬ 
ence electrode unless otherwise specified. 

The applicability of our electrochemical technique for measuring 
copper redistribution was examined in three different experiments. 
The first of these involved a cleaning and pretreatment procedure for 
AI 2024-T3. The Cu coverage was determined after each step of this 
procedure which included treatments in 0.1 M NaOH, 0.5 M 
HN0 3 + 0.65 M NaN0 3 , 3 M HN0 3 , and 0.1 M (NH 4 ) 2 S0 4 + 
0.03 M NH 4 OH. After each treatment the sample was washed with 
18 Mft Nanopure water and transferred to the three-electrode cell 
for elemental Cu coverage determination. In the second experiment, 
we examined the change in the amount of copper on the surface of 
an Al 2024-T3 sample as a function of time during immersion in 
0.1 M NaOH solution. In the third experiment, we examined the 
change in the amount of copper on the surface of an Al 2024-T3 
sample as a function of time during immersion in stagnant and 
stirred 0,5 M NaCl solutions. 

Results and Discussion 

Model experiment .—The model experiment described above was 
designed to determine the dependence of the associated Pb upd 
charge on the area of Cu in an Al-Cu sample. A relatively high elec¬ 
trolyte pH (pH 4) was chosen in order to suppress the hydrogen evo¬ 
lution reaction, that might be expected to occur in a high strength 
alloy, with a negligible overpotential. The millimolar amount of 
CH 3 COO~ present in the solution was expected to reduce the irre¬ 
versibility of the upd process and enhance the current peaks associ¬ 
ated with the Pb upd process. 15,16 

Figure 1 shows the current-voltage behavior of a polycrystalline 
copper electrode in solutions with (solid line) and without (dashed 
line) lead ions. The peak position (100 mV vs. Pb/Pb 2+ ) of the 
adsorption wave was close to that (105 mV) obtained by Schmidt 
et al. for the system Cu (polycrystalline)/ Pb 2+ in the presence of 
0.5 M KC1. The charge associated with the Pb upd layer was found 
to be 280 ± 10 jjlC cm -2 , about 65% of the one determined by 
Schmidt and Gygax. 11 This difference may be related to the coulom- 
bic contribution of Cl - ion adsorption on the copper surface, esti¬ 
mated to be 20-25% of the overall charge in the systems Cu (111)/ 



Figure 1 . Current-voltage behavior (cathodic sweep) of a Cu polycrystalline 
electrode in 10“ 1 M NaC10 4 + 10“ 3 M Pb(C10 4 ) 2 +10“ M HC10 4 + 
10“ 3 M CH 3 COONa solutions with and without 10“ 3 M Pb(C10 4 ) 2 . Sweep 
rate 20 mV s'" 1 . 


Pb 2+ , 10“ 4 Cl" 16 and Cu (001)/Pb 2+ , 10“ 4 Cl“. 17 These papers 
reported net charges of 300 and 262 jaC cm“ 2 for this system. 

Figure 2 shows the current-voltage behavior of Al-Cu samples 
containing different area fractions of Cu exposed in a solution con¬ 
taining Pb ions. These voltammograms were obtained using the fol¬ 
lowing background subtraction procedure. A polynomial form was 
fit to the dashed line (blank run containing no Pb ions) in Fig. 1 and 
this was subtracted from the solid curve of Fig. 1 yielding the 
Cu-100% wave shown in Fig. 2. The other curves shown in Fig. 2 
were obtained by numerically fitting the coefficients of the polyno¬ 
mial such that a separate "dotted line "or "blank run" for each sam¬ 
ple was generated that smoothly interpolated across the upd wave. 
These numerically generated blank runs have the advantage that a 
background subtraction can be accomplished without an a priori 
knowledge of the area fraction of the composite sample covered by 
copper. The results shown in Fig. 2 demonstrate that the peak height 
and area of the upd wave depend strongly on the area fraction of Cu 
exposed; the larger the area of Cu exposed the higher the amount of 
underpotentially deposited lead. 

Figure 3 is a normalized data collapse of the upd waves in Fig. 2 
showing that the wave shape and position is only very weakly depen¬ 
dent on the Cu content of the model Al-Cu sample. These compos¬ 
ite electrodes were polycrystalline in nature. For the samples in the 
range of 50-100% area fraction Cu, the upd lead sampled the identi¬ 
cal area of polycrystalline Cu, and within this composition range 
there was no shift in the upd peak position. For samples in the range 
of 0-50% Cu, the Cu portion of the composite surface was masked. 
Even so, we observed no significant (< —3 mV) shift in the peak 
position up to the sample containing 33% Cu. There was a shift, 
— 10 mV, for the samples containing 20 and 11% area fraction Cu. 
The polycrystalline Cu used in the background subtraction calibra¬ 
tion procedure had a grain size of order 200 fxm, so the 11% Cu sam¬ 
ple (that had an absolute Cu area of 0.11 cm 2 ) contained only about 
100 grains. The shift in the adsorption peak resulted from the upd 
process sampling a smaller number of grains for the samples con¬ 
taining less than —30% area fraction Cu. 

The final result extracted from the model experiment is present¬ 
ed in Fig. 4 which shows the normalized charges corresponding to 
the Pb underpotential wave as a function of the fractional area of Cu 
exposed. The best fit to the data (correlation coefficient = 0.998) is 
a straight line with a slope of 1. Comparison of the area fraction of 
Cu predicted by the linear fit to the actual values of the known stan¬ 
dards, showed that the amount of Cu exposed on the surface could 
be determined to an accuracy of ca. 2%. 



Figure 2. Current-voltage behavior of Al-Cu samples containing different 
area fractions of Cu exposed in 10“ 1 M NaC10 4 + 10“ 3 M Pb(C10 4 ) 2 + 
10“ 4 M HC10 4 + 10“ 3 M CH 3 COONa. Sweep rate -20 mV s -1 . 
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Figure 3. Normalized voltammetric curves for samples with different Al-Cu 
ratio in a solution containing 10" 1 M NaC10 4 + 10" 3 M Pb(C10 4 ) 2 + 
10" 4 M HC10 4 + 10~ 3 CH 3 COONa. Sweep rate -20 mV s" 1 . 


Pretreatment procedure for Al 2024-T3. —We examined the 
amount of copper remaining on the surface after each of four pre¬ 
treatment steps. The initial step was treatment of an aluminum alloy 
2024-T3 sample in 0.1 M NaOH solution at open-circuit potential 
for 1 min. This was expected to dissolve the gross oxides that may 
have formed on the alloy sheet surface during processing. Addition¬ 
ally this treatment should selectively leach Al from S-phase particles 
on the alloy surface. The polarizing step of the “copper removal” 
procedure proposed in Ref. 19 was used. The sample was treated 
accordingly in 0.5 M HN0 3 + 0.65 M NaN0 3 at a potential of 
-50 mV (SCE) for 40 min. The third step was a simple etching in 
3 M HN0 3 for 2 min at open-circuit potential, designed for the 
removal of any residual copper remaining on the surface. Finally, a 
procedure found to significantly reduce the cathodic activity of the 
alloy surface 20 was applied. The sample was treated in a solution 
containing 0.1 M (NH 4 ) 2 S0 4 + 0.03 M NH 3 (pH 8) at potentials 
-600 mV (SCE) and 150 mV (SCE) in order to clean the surface 
from all components that are expected to increase the surface 
cathodic activity, i.e., Cu, Fe, Mn. The results of this procedure are 
summarized in Fig. 5 where corrected voltammograms for all sam¬ 
ples examined are presented. The first curve presented is for lead 
upd on an untreated Al 2024 sample. After a treatment in a NaOH 
solution, a significant amount of Al had been selectively dissolved 



Figure 4. Normalized charges corresponding to the Pb underpotential wave 
as a function of the fractional area of Cu exposed. 



Figure 5. Corrected voltammetric curves obtained on Al 2024-T3 samples 
after treatment in different solutions (see text). Sweep rate -20 mV s" 1 . 


from the surface, resulting in a dramatic enrichment of elemental Cu 
on the alloy surface. The Pb upd result suggested that even after per¬ 
forming the copper removal procedure described in Ref. 19, there 
was still a significant amount of Cu remaining on the surface (Treat¬ 
ed with HN0 3 + NaN0 3 ). The most efficient treatment for copper 
removal appeared to be etching in 3 M HN0 3 for 2 min at open-cir¬ 
cuit potential. As a result of this treatment the residual copper on the 
surface (treated with HN0 3 :H 2 0) dropped below the originally 
existing level (Untreated Al 2024 T3). The final treatment did not 
show a significant decrease in the current signal coming from Pb upd 
[(NH 4 ) 2 S0 4 + NH 3 ]. 

Dealloying in NaOH solution .—Another illustration of the lead 
upd technique was performed examining the elemental copper con¬ 
tent of an Al 2024 T3 surface after treatment in 0.1 M NaOH at 
open-circuit potential for various times. The charge associated with 
the underpotentially deposited Pb layer after each treatment is plot¬ 
ted in Fig. 6. The results demonstrate that there is a dramatic 
increase of the elemental Cu on the alloy surface with increasing 
time. This is simply a result of progressive dealloying of Cu-con- 
taining intermetallics on the alloy surface. After the fifth minute of 
the experiment the charge exceeded a full monolayer level, implying 
that a significant “roughening” of a portion of the alloy surface 



Figure 6. Charge density of Pb upd layer formed on an Al 2024-T3 sample 
after treatment in 0.1 M NaOH solution for different times. 
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Figure 7. High-resolution scanning electron micrograph showing the mor¬ 
phology of the S-phase sample after polarization at 50 mV (SCE) in 0.1 M 
NaOH + A1C>2 solution. 
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Figure 8. Charge density of Pb upd layer formed on an A1 2024-T3 sample 
after treament in 0.5 M NaCl solution for different times. 


the surface of the alloy by a factor of about 1. 8 . The results discussed 
herein undoubtedly confirm the conjecture 1 that convection 
enhances the copper redistribution process. However, at present the 
precise mechanism of the enhanced redistribution associated with 
convection is unclear. 


occurred simultaneously with the copper enrichment process. This 
type of roughening is known to occur in alloy systems undergoing a 
dealloying process. During dealloying, the less noble component(s) 
of the alloy is selectively dissolved while the more noble element 
develops a bicontinuous solid-void morphology with a very large 
surface to volume ratio . 21 ' 23 We have recently shown that S-phase 
(Al 2 MgCu) intermetallic undergoes selective dissolution in NaOH 
and neutral chlorides resulting in a bicontinuous copper-void mor¬ 
phology . 24 A representative example for such a morpholgy is shown 
in Fig. 7. This high surface to volume area structure results in an ele¬ 
mental Cu surface coverage which can significantly exceed that of a 
single equivalent monolayer. Figure 6 shows that after 15 min of 
exposure there is a decrease in the charge of the Pb upd layer. We 
believe that this behavior is a result of ambient temperature coarsen¬ 
ing of the dealloyed morphology . 22 

Pitting in NaCl solution .—Using a procedure similar to that 
described above for NaOH, we examined Cu redistribution on the 
surface of A1 2024-T3 exposed to 0.5 M NaCl solution at open-cir¬ 
cuit potential for different times. It is well known that A1 2024-T3 
samples undergo pitting and intergranular corrosion 25,26 in a neutral 
chloride solution, and it was recently shown that that corrosion of 
this alloy in acid chloride results in a fourfold increase of the copper 
content on the surface of this alloy . 27 We believe that this enrichment 
in copper probably results from the selective dissolution of S-phase 
intermetallics. Additionally, if some of the nanometer-sized liga¬ 
ments of copper become mechanically disconnected from the bicon¬ 
tinuous dealloyed structure, these copper ligaments will dissolve at 
open circuit in the neutral chloride solution. The resulting copper 
cations will then plate back onto the alloy surface. Another explana¬ 
tion of the copper redistribution process associated with of A1 2024- 
T3 in a neutral NaCl solution was presented by Buchheit et al . 1 They 
suggested that small copper or copper-enriched remnants from the 
dealloying process mechanically detach from the surface and are 
mechanically dispersed and deposited onto the alloy surface by a 
non-Faradaic process, possibly involving convection. In our experi¬ 
ment, we sought to determine the effect of convection on the copper 
redistribution process. Figure 8 shows results for copper redistribu¬ 
tion obtained after treatments in stagnant and stirred (500 rpm) solu¬ 
tions. These results demonstrate that in these experiments convec¬ 
tion in the electrolyte increased the quantity of elemental copper on 


Summary 

A new surface-sensitive technique for determining the elemental 
copper coverage on the surface of Al alloys based on a upd process has 
been described. Model experiments on Al-Cu samples of known com¬ 
position demonstrated a linear dependence between the charge densi¬ 
ty associated with an underpotentially deposited Pb layer and the area 
fraction of Cu in the model system. An examination of different steps 
of a pretreatment procedures performed on Al 2024-T3 samples using 
this technique demonstrated that simple open-circuit etching of this 
alloy in 3 M HN0 3 to be most effective with respect to copper removal 
from the surface. Very significant enrichment of elemental Cu 
occurred when Al 2024 T-3 surfaces were exposed to 0.1 M NaOH 
solution or 0.5 M NaCl solution. Convection in the electrolyte was 
found to significantly enhance the Cu redistribution process. 
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Figure 1. Current-voltage behavior of Al, S phase, and Cu in 0.1 M NaOH + 
10~ 3 M AlOJ solution. Sweep rate -2 mV/s. 


the bulk of electrolyte) was estimated assuming D ~ 10 5 cm /s , 

5 ~ 5 x 10 " 2 cm , 9 n = 3, and F = 96,500 C mol -1 , to be 0.13 M 
dm" 3 . This value is in agreement with ancillary work that showed the 
solubility of NaA10 2 in a 0.1 M NaOH + 10 ~ 3 M A 1 2 (S 0 4 )3 elec¬ 
trolyte does not exceed 0.2 M dm" 3 . 

The polarization behavior of pure copper (Fig. 1) also shows an in¬ 
crease from a corrosion potential of 200 to '■'-TOO mV. At a po¬ 
tential of —100 mV we observed passivation-like behavior, followed 
by another increasing current region. The current voltage behavior is 
in general agreement with previously reported work of Strehblow and 
co-workers. For our experimental conditions (pH 13) this current 
wave can be attributed to CuO formation according to the reaction ’ 

lcu 7 0 + -H 2 0 - CuO + H + + e~ U1 

2 2 

The S-phase sample yields a dissolution curve with intermediate 
features of elemental Al and Cu. The open-circuit potential, 1.2 V, 
lies between that of aluminum and copper. From potentials of *1 
to -0.5 V, the current is almost potential independent. In this precrit- 
ical region, the alloy behaves as if it were “passive,” albeit with a rel¬ 
atively high current density. Figure 1 shows the critical potential 


behavior of the S phase, with the critical potential at ~ -450 mV. 
There is an additional current wave at ~ —100 mV, which corre¬ 
sponds to the wave observed on the copper electrode at +70 mV. 

Additional insight into the dealloying process can be obtained by 
analyzing the location of the peak in the current waves in the poten- 
tiodynamic curves for the S-phase sample and the pure copper. It is 
known that the electrochemical reactivity of the more reactive ele¬ 
ment in alloys is set by clusters of average size, £, which intersect the 
surface of the electrode . 12 The reactivity is expressed as the oxida¬ 
tion potential defined by 

AE = E{p) - E = ^ [2] 

nq k 

where 7 is the liquid/solid interfacial energy, ft is the atomic volume, 
n is the number of electrons transferred in the dissolution process, and 
q is the elementary charge of an electron. E is the reversible potential 
of flat elemental A defined by Nemst equation. Using the relation 
above, we can estimate the average initial size of the copper liga¬ 
ments resulting from S-phase dealloying. For 7—1 J/m 2 , ft - 5 X 
10-29 m 3 } q = 16 x _KT 19 C, n = 1 (Eq. 1) and measuring the 
potential shift [E(p) - E] = 0.170 V from Fig. 1, the average initial 
(prior to coarsening) metal ligament diameter is of order 2 nm. 

The curve in Fig. 3 shows the current-time behavior of the S phase 
sample at a potential of -200 mV where Al is being selectively dis¬ 
solved. The current shows oscillatory behavior, which has been ob¬ 
served in other systems and is related to a complex interaction be¬ 
tween metal dissolution and mass transport effects in the electrolyte. 
As aluminum dissolves from the S phase above the critical potential, 
a porous structure develops with pits of order 2 nm in diameter. Metal 
dissolution within these pits is inhibited by mass transport. The oscil¬ 
lations are related to the stability of a film (possibly a salt film, 
NaA10 9 , or a hydrated oxide) at the bottom of the pits. Film stability 
is determined by a delicate balance between metal dissolution, ion 
solubility, and diffusion transport in the electrolyte. As dissolution 
proceeds, the electrolyte becomes saturated in the aluminate anion, 
which results in a precipitation event causing the drop in the current 
owing to film resistance. Mass transport in the electrolyte results in a 
decrease in the concentration of A10 2 to below saturation at the salt 
film/electrolyte interface. This allows for the chemical dissolution of 
the film. At this stage the metal starts to dissolve again until the elec¬ 
trolyte near the surface of the metal within the pits is once again sat¬ 
urated in the metal anion and precipitation occurs again. This process 
occurs in parallel for every pit in the porous structure. Since tine trans¬ 
port fields of the pits are coupled, complex nonlinearities develop. 
Such systems are known to display chaotic behavior . 13 



.2 1-1-1-1-1- 1 - 1 - - - 1 - 1 - 
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Figure 2. Curve fit of the potentiodynamic curve for the system Al/0.1 M 
NaOH + 10 -3 M AIO 2 with the function ti = const (1 - ilif). Correlation 
coefficient, 0.99. i 4 = 7.92 mA/cm 2 . 



Figure 3. Current-time curve during long-time corrosion of an S-phase sam¬ 
ple at a potential of —200 mV. 
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Figure 4. SEM showing the structure of an S-phase sample after 20 h of 
immersion in 0.1 M NaOH 4- 10~ 3 M AIO 2 solution at a potential of 
-800 mV, which is below the critical potential. 



Figure 6. High-magnification SEM showing the structure of S-phase sample 
after 20 h of immersion in 0.1 M NaOH + 10“ 3 M AlOJ solution at a poten¬ 
tial of —200 mV. 


High resolution scanning electron microscopy (SEM).— Follow¬ 
ing dealloying of the S phase samples we used high resolution SEM 
to examine the morphology of the surface. Figure 4 shows the mor¬ 
phology of the surface obtained after 20 h of immersion in 0.1 M 
NaOH + 10~ 3 M AIO 2 solution at a potential of —800 mV (below 
the critical potential). The micrograph reveals some morphology 
evolution related to corrosion, but the bicontinuous microstructure 
characteristic of a bulk dealloying process 7 does not evolve. Figures 
5 and 6 show a bicontinuous dealloying morphology, that resulted 
from the selective dissolution of S phase after 20 h of immersion in 
0.1 NaOH + 10~ 3 M A10 2 ” solution at a potential of —200 mV 
(above the critical potential). The flatter regions showing “Swiss 
cheese’Mike morphology are typical of room temperature coarsen¬ 
ing and sintering processes in dealloying microstructures. 7 Sintering 



Figure 5. SEM showing the structure of an S phase sample after 20 h of 
immersion in 0.1 M NaOH + 10~ 3 M AIO 2 solution at a potential of 
-200 mV which is above the critical potential. 


results from mass-induced collapse of the nanoscale copper-rich lig¬ 
aments in the near surface region of the porous structure. Beneath 
this surface layer there is a bicontinuous solid-void microstructure 
with length scale of order 20-50 nm. This microstructure is reminis¬ 
cent of the appearance of Cu-Au and Ag-Au surfaces, which have 
undergone dealloying under appropriate conditions . 7 

Conclusions 

We present electrochemical and SEM evidence of S-phase deal¬ 
loying in an alkaline solution. The current-voltage curve obtained 
demonstrates critical potential behavior and selective dissolution un¬ 
der the experimental conditions examined. The current-time curve 
obtained after long-term dealloying suggested a complex interaction 
between metal dissolution and mass-transport effects in the electro¬ 
lyte, which resulted in what is probably chaotic behavior of the sys¬ 
tem investigated. S-phase corrosion above the critical potential re¬ 
sulted in the formation of bicontinuous microstructures, typical of 
dealloying processes. 
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mass transport through the bulk of the growing phase (the electro- 
lyte) is always a stabilizing influence 12 , and because mass transport 
through the bulk of the dissolving phase appears too slow to be 
significant, we hypothesized that the morphology-determining 
physical process is confined to the interface region between the 
alloy and the electrolyte. To test this, we developed a kinetic Monte 
Carlo model to simulate Ag-Au dealloying, including only diffusion 
of silver and gold and dissolution of silver 13 . We found that this 
model was able to reproduce all relevant experimental trends 
characteristic of dealloying, both morphological and kinetic. 

Figure 2 shows a simulated porous structure with ligament widths 
of 2-5 nm. Our simulations were successful in modelling the 
nanoporous morphology, and also in modelling the dynamic 
behaviour of the dissolution current versus overpotential. It is a 
well characterized feature of alloy dissolution that as the over¬ 
potential is increased (usually at rates of the order of a few millivolts 
per second), the dissolution current of ions from the alloy stays at a 
low level until a bulk-composition-dependent critical potential (V c ) 
is reached, at which point this current rises rapidly 14 . Figure 3 
shows simulated and experimental polarization curves for different 
alloy compositions. There is clear observation of a composition- 
dependent V c . To our knowledge, this is the first simulation model 
to produce such behaviour, suggesting that we have found a 
minimum set of physical processes to include in any model for 



Figure 2 Simulated nanoporous gold. The simulation model was as follows: a bond¬ 
breaking model was used for diffusion; atoms with N near neighbours diffused with rate 
k N = v D exp(—A/e/Ag7), where e is a bond energy and v D = 10 13 s“ 1 . Dissolution rates were 
consistent with the Butler-Volmer (BV) equation in the high-driving-force Tafel regime; 
the dissolution rate ^for a silver atom with N near neighbours was written as k EiN = 
v £ expH/Ve - where v E = 10 4 s" 1 is an attempt frequency determined by the 
exchange-current density in the BV equation and <f> is the overpotential. For the figure, 
0 = 1.75 eV, elk^T- 5.51. 


alloy dissolution. 

The simulations reveal the following qualitative picture of 
porosity formation. The process starts with the dissolution of a 
single silver atom on a flat alloy surface of close-packed (111) 
orientation, leaving behind a terrace vacancy. The atoms coordinat¬ 
ing this vacancy have fewer lateral near-neighbours than other silver 
atoms in the terrace, and are thus more susceptible to dissolution. 
As a result, the entire terrace is ‘stripped’, leaving behind gold atoms 
with no lateral coordination (‘adatoms’). Before the next layer is 
attacked, these gold adatoms diffuse about and start to agglomerate 
into islands. As a result, rather than a uniform diffuse layer of gold 
spread over the surface, the surface is comprised of two distinct 
kinds of regions—pure gold clusters that locally passivate the 
surface, and patches of un-dealloyed material exposed to electrolyte. 
When silver atoms in these patches dissolve, more gold adatoms are 
released onto the surface. These adatoms diffuse to the gold clusters 
left over from dissolution of previous layers, continuing to leave un- 
dealloyed material exposed to electrolyte. In the early stages, these 
gold clusters are mounds that are gold-rich at their peaks and that 
have alloy composition at their bases. These mounds get undercut, 
increasing the surface area that gold must cover to bring about 
passivation. Ultimately, this leads to pit formation and porosity. 

Central to this description is the the coalescence of gold adatoms 
into stable clusters. The spacing between these ‘islands’ in the initial 
stages of dissolution is close to the spacing between ligaments in the 
final porous structure. The physical reason for this coalescence can 
be understood by considering the gold adatoms to be one compo¬ 
nent of a two-component solution of gold and electrolyte confined 

to the monolayer-thick interfacial layer sitting on top of un- 
dealloyed material. We modelled the thermodynamics of the inter¬ 
facial layer as a regular solution 15 , and found the solubility of gold in 
electrolyte within the interfacial layer to be of the order of 10 per 
site (see Methods). This solubility may be interpreted as the 
‘equilibrium concentration of gold adatoms on the surface of the 
alloy—in the absence of etching, it represents a dynamic equilib¬ 
rium of adatoms resulting from their two-dimensional evaporation 
from step edges onto terraces and their subsequent recondensation. 

In contrast to the equilibrium condition, rapidly stripping a 
terrace of silver atoms leaves gold adatoms with a local site 
occupancy fraction equal to that in the bulk, typically 10-40%— 
far above their equilibrated concentration of 10 7 per site. Thus, 
there is an extremely strong driving force for gold adatoms to 
condense onto nearby gold-rich clusters. In fact, regions of surface 
with high enough supersaturation of gold adatoms sit ‘within the 
spinodal’, a special segment of the curve representing free energy / 
of a spatially uniform layer versus gold concentration C for which 
d 2 f IdC 2 < 0. Within the spinodal, composition fluctuations of 
infinitesimal amplitude lead to a lower overall free energy for the 




Figure 3 Comparison of experimental and simulated current-potential behaviour, 
a, Current-potential behaviour tor varying Ag-Au alloy compositions (atom% Au) 
dealloyed in 0.1 M HCI0 4 + 0.1 M Ag + (reference electrode 0.1 M Ag + /Ag). b, Simulated 


current— potential behaviour of Ag-Au alloys, c, Comparison of experimental (line) and 
simulated (triangles) critical potentials; the zero of overpotential has been set equal to the 
onset of dissolution of pure silver both in simulation and in experiment. 
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system, and involve atomic diffusion against concentration gradi¬ 
ents (the uphill diffusion’ process through which gold condenses 
onto nearby clusters), that is, the system is inherently unstable and 
will spontaneously phase-separate. But fluctuations of long length 
scale grow slowly due to the required diffusion times, and short- 
length-scale fluctuations create much energetically unfavourable 
incipient interface between the phases, inhibiting their growth. 
Hence, phase separation is manifested most rapidly at an inter¬ 
mediate length scale that roughly corresponds to the spacing 
between the observed gold-rich clusters. This effect is known as 
spinodal decomposition 16,17 . As porosity forms, the decomposition 
is occurring on a non-flat, non-uniform surface with continuously 
increasing surface area. 

In our model, the motion of the alloy-electrolyte interface is fully 
described mathematically by the flux of diffusing adatoms / s , the 
velocity of the interface normal to itself v n , and the concentration 
accumulation rate dCldt , all of which are interrelated and vary with 
position along the curve of the interface (for detailed derivations, 
see Methods). For J s , we used a model for diffusion during spinodal 
decomposition known as the Cahn-Hilliard equation 17 . The 
normal velocity depends on C and also on the local curvature k 
through capillary effects 11 . The time evolution of C is uniquely 
determined by the local mass-conservation condition 

dC/dt - v n C 0 ~ v n KC - V-/ s (1) 

where C 0 is the bulk gold concentration. This condition is analogous 
to the local conservation of heat or solute appearing in boundary- 
layer models of solidification 18 , with two important differences: 
( 1 ) the interfacial layer thickness is constant along the interface, and 
is microscopic, rather than being a spatially varying macroscopic 
diffusion length, and ( 2 ) the surface aggregation process inherent in 
the Cahn-Hilliard form for J s is essential for porosity formation. 
Simple (‘downhill’) surface diffusion (J$ = —D$VC) yields an 



Figure 4 Simulated evolution of an artificial pit in Au 10 <%Ag 90 % (atom%), <£ = 1.8 eV. 
Cross-sections along the (11 i) plane defined by the yellow line in a are shown below each 
plan view, a, The initial condition is a surface fully passivated with gold except within a 
circular region (the 'artificial pit’), b, After 1 s, the pit has penetrated a few monolayers into 
the bulk, We note how there are fewer gold clusters near the side wall than at the centre of 
the pit. c, After 10 s, a gold cluster has nucleated in the centre of the pit. d, At 100 s, the 
pit has split into multiple pits; each will continue to propagate into the bulk to form a 
porous structure like that in Fig. 2. 


initially unstable interface that passivates quickly, before well- 
formed pores have a chance to develop. Cahn-Hilliard diffusion 
also dominates capillarity-driven surface diffusion—the effect 
usually incorporated into interface evolution equations 19 . 

We performed numerical integration of equation (1) using a 
relative arc-length parametrization scheme 20,21 , and parameters that 
matched those used in the kinetic Monte Carlo simulations. We 
observed, as expected, the evolution of gold clusters separated by a 
characteristic spacing X. An analytic expression for X can be found 
by a time-dependent linear stability analysis of equation (1) that 
takes into account the slow increase of gold concentration into the 
interfacial layer as the instability develops. This effect needs to be 
included because the spatial period with the largest amplification 
rate depends sensitively on the gold concentration. Specifically, this 
spatial period decreases sharply as the concentration increases past a 
threshold concentration for instability that corresponds to the 
spinodal point d 2 f IdC 2 = 0, and the interface is stable for concen¬ 
trations below this threshold. We find a maximally unstable spatial 
period that scales as X ^ (D s /Vo) 1/6 > where V 0 is the velocity of a flat 
alloy surface with no gold accumulated upon it. This prediction is in 
qualitative agreement with both kinetic Monte Carlo simulations 
and experiments, both of which show that the characteristic length 
scale of porosity decreases with increasing driving force. A more 
elaborate analysis incorporating nonlinear effects, however, remains 
needed for a detailed quantitative comparison. 

There is an analogy between this result, applicable to etching, and 
two-dimensional island nucleation during submonolayer vapour 
phase deposition. Namely, in the early stages of etching, the 
dissolution process is analogous to deposition of gold; in both 
processes, adatoms are added to the surface where they are free to 
agglomerate into islands. The case of vapour phase deposition has 
been studied using rate equations that describe an aggregation 
process where adatoms stick together irreversibly. In these studies 
it is a well-known result that the island spacing scales as (D s /jF)' i , 
where the deposition rate F is the direct analogue of the surface 
velocity in etching, and the exponent (x depends on details of the 
aggregation process (see ref. 22 and references therein). That these 
results are limited to irreversible aggregation during deposition and 
our analysis is for reversible aggregation during etching suggests the 
existence of universal scaling laws for aggregation that do not 
depend on reversibility or the lack thereof in these two opposite 
processes. 

Our kinetic Monte Carlo simulation elucidates the later stages of 
morphological evolution, and the mechanism by which three- 
dimensional porosity evolves. We highlight the features of this 
process by showing in Fig. 4 a simulation of an artificial pit in an 
otherwise fully passivated surface. When the pit reaches sufficient 
depth, its surface area has increased sufficiently that a new gold 
cluster nucleates. When this happens, the pore splits into multiple 
new pits, each with a smaller surface area than its parent. These 
‘child’ pits continue to penetrate into the bulk, increasing their 
surface area, nucleating new clusters, spawning new pits, and so on, 
until a full, three-dimensional porous structure evolves, such as 
those illustrated in Figs 1 and 2. □ 

Methods 

In a regular solution, the enthalpy of mixing depends on the bond energies and the entropy 
of mixing is ideal. The free energy of a regular solution f(C,T) is written / = ac(l - c)+ 
fc B T[cInc + (1 - c)ln(l “ c)], where c is the mole fraction of gold (c = CO 2 *, where O is 
atomic volume), a = 6[£ Au -rieclrnlylc — (1/2 )(£a«-Au "l" •£eleclr«fyte-elcciro!yle)]> where H x -x the 
respective interaction energies between Au and electrolyte, the prefactor 6 is the lateral 
coordination in the two-dimensional hexagonal lattice of the interfacial layer, is 
Boltzmann’s constant and T is absolute temperature. For our simulation conditions, 
realistic timescales and length scales were obtained from the parameters 
£ Au _ Au = -0.285 eV (= -e, the simulation bond energy as described in Fig. 2), T = 600 K, 
^Au-electrolyte = Edcoroiyr-dcctmiytc = 0-0 eV. With these parameters, a = 0.855 eV. The free 
energy has the familiar double-well form” and a minimum at c ~ 10” per site, 
representing the solubility of gold in electrolyte (and vice versa). 

The Cahn-Hilliard diffusion equation is 7 S = -M{C)(d 2 fJdf)VC + 2M(QwV , C. 
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Here, 453Af(C) is a mobility, w is the so-called gradient energy coefficient, and the 
gradients are taken with respect to arc length. The first term on the right-hand side 
describes the chemical effect leading to phase separation within the spinodal; the second 
term describes the effect that damps short-wavelength fluctuations. The mobility is 
proportional to the surface diffusivity D s and is given by M(C) = {Ds/k^T)c(l-c). The 
mobility is peaked for c - 0.5, and is zero for c - 0 and c = 1 (atoms do not diffuse in pure 
phases because there are no vacancies in our model). The normal velocity is given by 
v„(C) — V(C)[] - ( yQlk H T )#c], where 7 is the surface free energy and V(C) is called the 
interface response function, equal to the velocity of a flat surface covered with a 
concentration C of gold. We find in both simulation and experiment that the interface 
response is fitted well by the functional form V(C) = V 0 (^)exp(-C/C*), where 4> is the 
overpotential and C* is a constant. Experimentally, the gold accumulation can be inferred 
by integrating the dissolution current versus time at fixed overpotentiai; it is necessary to 
use an overpotentiai that is low enough to ensure that the surface remains planar (that is, 
porosity does not form) and also to catch the short initial transient rise in current as silver 
atoms are pulled from the first few monolayers. This particular form for the interface 
response function is quite curious. Naively, one might expect that the local interface 
velocity would be proportional to the local concentration of silver exposed to the 
electrolyte, that is, V(Q « (1 -c). However, the decaying exponential form suggests that 
there is an evolving distribution of holes opening and closing within the interfacial region, 
controlling the accumulation rate. 

Physically, the mass conservation condition (equation (1)) is the statement that the 
total number CbAs of gold atoms in a length As of interface with lateral width b can 
change as a result of three distinct effects that correspond to the three terms on the 
right-hand-side of equation ( 1 ): the accumulation of gold atoms into the interfacial 
layer from the solid being dissolved; the local stretching of the interface (cJAs/df = v„k 
A s), which can either increase or decrease C depending on whether the solid is concave 
(k > 0 ) or convex (<c < 0 ); and the motion of atoms along the interface driven by the 
surface diffusion flux J s . 

Received 14 November 2000; accepted 10 January 2001. 

]. Pickering, H. W. Characteristic features of alloy polarization curves. Corros. Sci. 23, ) 107—1120 
(1983). 

2. Forty, A. J. Corrosion micromorphology of noble metal alloys and depletion gilding. Nature 282, 597- 
598 (1979). 

3. Masing, G. Zur Theorie der Resistenzgrenzen in Mischkristallen. Z. Anorg. AUg. Chem. 118, 293-308 
(1921). 

4. Lechtman, H. Pre-Columbian surface metallurgy. Sci. Am. 250, 56-63 (1984). 

5. Williams, D. E., Newman, R. C., Song, Q. 8c Kelly, R. G. Passivity breakdown and pitting corrosion of 
binary alloys. Nature 350, 216-219 (1991). 

6. Newman, R. C. & Sieradzki, K. Metallic Corrosion. Science 263, 1708-1709 (1994). 

7. Li, R. & Sieradzki, K. Ductile-brittle transition in random porous Au. Phys. Rev. Lett. 68, 1168-1171 
(1992). 

8. Corcoran, S. in Critical Factors in Localized Corrosion III (eds Kelly, R. G., Frankel, G. S., Natishan, 

P. M. & Newman, R. C.) 500-507 (Electrochemical Society, Pennington, New Jersey, 2000). 

9. Pickering, H. W. & Wagner, C. Electrolytic dissolution of binary alloys containing a noble metal. 

/. Elcctrochem. Soc. 114, 698-706 (1967). 

10. Sieradzki, K., Corderman, R. R., Shukla, K. 8c Newman, R. C. Computer simulations of corrosion: 
selective dissolution of binary alloys. P/ii7. Mag. A 59, 713-746 (1989). 

11. Sieradzki, K. Curvature effects in alloy dissolution. }. Elcctrochem. Soc. 140, 2868-2872 (1993). 

12. Wagner, C. Contribution to the theory of electropolishing. J. Bectrochem. Soc. 101, 225-228 
(1953). 

13. Erlebachcr, J. in Dynamics of Crystal Surfaces and Interfaces (eds Duxbury, P. & Pence, T.) 24-28 
(Plenum, New York, 1997). 

14. Wagner, K., Brankovic, S. R., Dmitrov, N. 8c Sieradzki, K. Dealloying below the critical potential. 

}. Elcctrochem. Soc. 144, 3545-3555 (1997). 

15. Cahn, J. W. & Hilliard, J. E. Free energy of a nonuniform system. 1. Interfacial free energy. /. Chcm. 
Phys. 28, 258-267 (1958). 

16. Cahn, J. W. & Hilliard, J. E. Free energy of a nonuniform system. III. Nucleation in a two-component 
incompressible fluid. /. Chem. Phys. 31, 688-699 (1958). 

17. Hilliard, J. E. in Solidification 497-560 (American Society for Metals, Metals Park, Ohio, 1971). 

18. Ben-Jacob, E„ Goldenfeld, N., Langer, J. S. 8c Schon, G. Dynamics of interfacial pattern formation. 
Phys. Rev. Lett. 51, 1930-1932 (1983). 

19. Mullins, W. W. Theory of thermal grooving. ]. Appl. Phys. 28, 333-339 (1957), 

20. Kessler, D. A., Koplik, J. & Levine, H. Geometrical models of interface evolution. II. Numerical 
simulation Phys. Rev. A 30, 3161-3174 (1984). 

21. Hou, T. Y„ Lowengrub, J. S. 8c Shelley, M. J. Removing the stiffness from interfacial flow with surface 
tension. ]. Comp. Phys. 114, 312-338 (1994). 

22. Barabasi, A. -L. 8c Stanley, H. E. Fractal Concepts in Surface Growth (Cambridge Univ. Press, 1995). 

23. Haasen, P. Physical Metallurgy, (Cambridge Univ. Press, 1986). 

24. Forty, A. J. in Sir Charles Frank An 80th Birthday Tribute (ed. Chamber, R .G.) 164-187 (Adam Hilger, 
Bristol, 1991). 

25. Tulimieri,D. J., Yoon, J. 8c Chan, M. H.W. Ordering of helium mixtures in porous gold. Phys. Rev. Lett. 
82, 121-124 (1999). 


Acknowledgements 

This work was supported by the US Department of Energy, Basic Energy Sciences. 
The research of A.K. also benefited from computer time allocation at NU-ASCC. 
K.S. thanks the AFOSR for support. 

Correspondence and requests for materials should be addressed to J.E 
(e-mail: Jonah.Erlebacher@jhu.edu). 


Ice shelves in the Pleistocene Arctic 
Ocean inferred from glaciogenic 
deep-sea bedforms 

Leonid Polyak*, Margo H. Edwardst, Bernard J. Coakley* 

& Martin Jakobsson§!! 

* Byrd Polar Research Center , Ohio State University, Columbus, Ohio 43210, USA 
t Hawaii Mapping Research Group, Hawaii Institute of Geophysics and 
Planetology, University of Hawaii, Honolulu, Hawaii 96822, USA 
$ Department of Geology, Tulane University, New Orleans, Louisiana 70118, USA 
§ Department of Geology and Geochemistry, Stockholm University, 

106 91 Stockholm, Sweden 


It has been proposed that during Pleistocene glaciations, an ice 
cap of 1 kilometre or greater thickness covered the Arctic 
Ocean 1 ’ 3 , This notion contrasts with the prevailing view that the 
Arctic Ocean was covered only by perennial sea ice with scattered 
icebergs 4 " 6 . Detailed mapping of the ocean floor is the best means 
to resolve this issue. Although sea-floor imagery has been used to 
reconstruct the glacial history of the Antarctic shelf 7 " 9 , little data 
have been collected in the Arctic Ocean because of operational 
constraints 1011 . The use of a geophysical mapping system during 
the submarine SCICEX expedition in 1999 12 provided the oppor¬ 
tunity to perform such an investigation over a large portion of the 
Arctic Ocean. Here we analyse backscatter images and sub-bottom 
profiler records obtained during this expedition from depths as 
great as 1 kilometre. These records show multiple bedforms 
indicative of glacial scouring and moulding of sea floor, combined 
with large-scale erosion of submarine ridge crests. These distinct 
glaciogenic features demonstrate that immense, Antarctic-type 
ice shelves up to 1 kilometre thick and hundreds of kilometres 
long existed in the Arctic Ocean during Pleistocene glaciations. 

The central Arctic Ocean contains relatively shallow areas (water 
depths < 1,000 m; see Fig. 1) on Yermak plateau, Lomonosov ridge 
and Chukchi borderland—which includes Chukchi plateau, Chuk¬ 
chi rise and Northwind ridge. During the SCICEX-99 expedition, 
conducted on the nuclear-powered submarine USS Hawkbill , 
shallow sea-floor areas were targeted for mapping to detect glacio¬ 
genic bedforms. Sea-floor images (collected using a submarine- 
mounted 12-kHz swath bathymetry and sidescan sonar 12 ) from the 
Chukchi borderland and the Lomonosov ridge show a variety of 
bedforms, including random or subparallel scours, parallel linea- 
tions, and transverse ridges. On the records from the chirp sub¬ 
bottom profiler, these bedforms are associated with planed ridge 
crests with rough microrelief and obvious angular unconformities 
cut into the stratified sediments. 

Randomly oriented furrows, typically <100-m wide and up to 
30-m deep, densely cover the shallowest, <400-m-deep portions of 
sea floor on the Chukchi borderland and adjacent continental 
margin (Fig. 2a). Isolated larger scours up to 700-m wide and 
over 10-km long occur as deep as 500 m. Even greater depths, 
exceeding 900 m, are attained by closely spaced, subparallel scours 
on the Lomonosov ridge. Sea-floor scours are known to be formed 
by the drift of icebergs and pack-ice ridges 13 . At present, icebergs in 
the Arctic Ocean have at most 50-m draughts 14 , whereas icebergs off 
Antarctica and Greenland reach depths of 500-550 m (refs 15, 16). 
The largest depths of gouged sea floor, extending to 850 m, have 
been reported from the Yermak plateau 10 , matching the depth of 
sours on the Lomonosov ridge. 

Below the depth range of dense scouring, the sea floor exhibits 
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We present results of a comprehensive electrochemical study aimed at determining the factors that control the corrosion evolution 
of the important structural alloy A1 2024-T3. An important aspect of corrosion in this alloy is copper redistribution. The redis¬ 
tributed copper serves to enhance the kinetics of the corrosion process and presents some major difficulties with respect to 
development of new families of conversion coatings. Two major sources of “mobile” copper are discussed, the Al 2 CuMg S-phase 
and copper in solid solution. In each case, a dealloying mechanism is responsible for liberating copper that becomes available to 
redistribution by solid state and/or liquid phase transport mechanisms. This work is concerned with identifying the relative 
importance of these sources of copper and the active transport mechanisms. In order to decouple the sources of copper, we 
designed "synthetic” 2024 alloys that emulate the oxygen reduction kinetics of the real alloy. Our results indicate that, in 
principal, both sources of copper are operative during redistribution and, for the particular alloy that we examined, both contribute 
about equally to the amount of redistributed copper. 

© 2002 The Electrochemical Society. [DOI: 10.1149/1.1498258] All rights reserved. 
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Copper redistribution occurs on the surface of copper containing 
aluminum alloys such as 2024-T3 as a result of dealloying. There 
are two principal sources of copper in the alloy: intermetallic copper 
containing precipitate particles such as the orthorhombic S-phase 
(Al 2 CuMg) in 2024-T3, and copper in solid solution. 1 ' 4 Henceforth, 
we refer to dealloying of matrix aluminum and the buildup of cop¬ 
per from solid solution as matrix dealloying. The surface coverage 
of S-phase particles on 2024-T3 is of order 0.03 area fraction, 1 and 
there is —0.2 to 0.5 wt % Cu in solid solution. 5 Assuming a mean 
S-phase particle diam of 1 jxm, the mean interparticle separation is 
—5 fxm. 

Figure 1 is a schematic illustration of S-phase dealloying which 
occurs in corrosive environments containing, for example, chlorides, 
or a commercial deoxidizing pretreatment solution. Various re¬ 
searchers have also noted similar effects during polishing of these 
alloys in nominally pure water. 6 ' 9 The oxide on the surface of the 
S-phase is not very protective and a fluctuation (transient break¬ 
down event) or heterogeneity in the oxide exposes bare metal to the 
electrolyte. Once this event occurs, Mg dissolves into the electrolyte 
as the soluble Mg 2+ cation. The dynamics of this process transiently 
inhibits repassivation and this allows for aluminum dissolution from 
the S-phase. The remaining Cu forms a very high surface area po¬ 
rous sponge structure as shown in Fig. 2, which is typical of bicon- 
tinuous dealloyed microstructures. 10 Thus the S-phase is initially 
anodic, but after dealloying, turns into a high surface area cathode 
for oxygen reduction. As long as the remnant Cu sponge remains 
mechanically and electrically connected to the A1 matrix the corro¬ 
sion potential is too low to allow for Cu dissolution by electro¬ 
chemical means. However, if the particles of the Cu sponge are 
mechanically disconnected (Fig. lb-d), the copper particle will dis¬ 
solve at its corrosion potential in, for example, a Cl - electrolyte, and 
electrochemically plate back onto the aluminum alloy surface. We 
have performed separate tests which confirm this scenario. The elec¬ 
trodeposition of dissolved copper onto the oxide covered alloy sur¬ 
face may at first seem surprising, but there are regions of this sur¬ 
face that are at least transiently bare. Additionally, the oxide 
on the alloy surface is undoubtedly highly defective and a 
semiconductor. 11 

In a neutral chloride solution, oxygen reduction occurs according 
to 2H 2 0 + 0 2 + 4e“ = 4(OH)“, which results in an increase in 
the pH of the electrolyte in the neighborhood of local cathodes. The 
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Cu-rich remnant S-phase sponge serves as cathode for this reaction. 
As the pH increases to about the range 9-9.5, 8,12,13 the passive oxide 
on the surface of the alloy matrix will chemically dissolve and the 
bare aluminum will begin to selectively dissolve via the soluble 
AlOj anion. This scenario is illustrated in Fig. 3. As hydroxide is 
produced at a local cathode, it diffuses radially outward activating 
an ever increasing surface area surrounding the cathode to alkaline 
dissolution of aluminum. This matrix dealloying process results in 
the agglomeration of solid solution Cu into Cu particles that also 
serve as cathodic sites available to support yet more oxygen reduc¬ 
tion. Additionally, dissolution of aluminum will bring new precipi¬ 
tate particles to the surface, which can activate anodically and sub¬ 
sequently turn in to reduction sites. This situation is complicated 
further, since the cathodes have finite lifetimes. A copper particle 
can stop functioning as a cathode either because it eventually is 
undercut by alkaline dissolution of aluminum and becomes discon¬ 
nected, or from buildup of a gel-like aluminum hydroxide which 
blocks the cathode from access to oxygen. 

Figure 4 is a schematic illustration of the growth of the zone of 
increased pH in the vicinity of a dealloyed S-phase particle. In the 
long-term limit, the growth of the zone will be dictated by the ge¬ 
ometry of overall mass transport in the electrolyte, either radial dif¬ 
fusion of 0 2 toward the cathode or (OH)" diffusion radially away 
from the cathode. Each dealloyed particle behaves as a microelec¬ 
trode, since the particle diameter is considerably smaller than the 
width of the diffusion boundary layer. Thus, the zone of increased 
pH is initially planar (when the zone size is of the order of the 
cathode diameter) with Cottrellian current but eventually adopts the 
hemispherical growth front that characterizes the diffusive field. The 
zone grows in this geometry until its radius is of the order of the 
diffusion boundary layer thickness (10-100 p,m) at which point the 
growth of the zone is cut off. Thus the diffusion boundary layer 
thickness sets an upper limit on the radius of the increased pH zone 
that can be developed by a single particle. 

If the mean S-phase particle diameter is 1 fim, then the average 
separation between particles will be ~5 \xm corresponding to an 
—0.03 area fraction surface coverage of copper. Figure 5 illustrates 
the growth and overlap of these zones of increased pH emanating 
from an array of these cathode microelectrodes. If the mean particle 
spacing is less than the diffusion boundary layer thickness, the hemi¬ 
spheres of increased pH overlap. As this occurs, the diffusion field 
switches from radial to a one-dimensional (1-D) field in the direc¬ 
tion of the surface normal and the current becomes Cottrellian. In 
this case, the entire surface is activated to matrix dealloying and a 
significant increase in the surface copper concentration develops. 
This situation is physically similar to that of a nucleation and growth 
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Figure 1. Schematic illustration showing dealloying of an S-phase particle, (a) The quenched-in atomic structure of an S-phase particle prior to dealloying. 
Only the atomic-scale copper in this structure is shown, (b) As the aggressive fluid phase removes Mg and Al, the remnant Cu forms a nanometer-scale porous 
sponge structure part of which is mechanically disconnected from the solid surface in the electrolyte, (c) As the porous Cu structure evolves more Cu particles, 
on the order of tens of nanometers in diameter, are liberated and (d) is a continuation of this process. Box size length scale, 1 jxm. 


process with overlapping growth fronts. It is known that such a 
process is described by Avrami growth kinetics; thus the partial 
current of either oxygen reduction or aluminum dissolution should 
be given by an equation of the form 7=1- exp(-af”). Here, 7 is 
the normalized current density (normalized to the value characteris¬ 
tic of a fully activated surface), a is a constant related to the geom¬ 
etry of the extending zone of increased pH, and n corresponds to 
1 + the dimension of the growth assuming a continuous nucleation 
process. This situation is, however, complicated by the fact that 
cathodes are dying out, and the “nucleation time'’ of the cathode is 
controlled by fluctuations. Thus a proper formulation of this prob¬ 
lem must include terms that correspond to both cathode nucleation 
processes and cathode extinction processes. 

This paper is aimed at addressing the relative contributions of 
each of the mechanisms described above to copper redistribution on 
the surface of Al alloy 2024-T3. The manuscript is organized as 
follows. In the next section we describe the general experimental 
methods used in this study. Then we present results from rotating 
disk electrode (RDE) and rotating ring-disk electrode (RRDE) stud¬ 
ies for copper redistribution on the surface of Al alloy 2024-T3 (real 



Figure 2. SEM showing the structure of a dealloyed S-phase sample after 20 
h of immersion in 0.1 M NaOH + 10 -3 M AlO^ solution. 


2024) and synthetic 2024-T3 (synthetic 2024) alloys. Finally, we 
present results for oxygen diffusion-limited current densities on 
RDEs and in the conclusions, we knit together all the results in a 
model that presents a coherent qualitative picture of corrosion and 
copper redistribution on the surface of Al alloy 2024-T3. 

Experimental 

Sample preparation. —We used five different kinds of samples in 
our study, Al 2024-T3 alloy, Al-1 wt % Cu solid solution alloys, 
synthetic 2024 (see the section on Synthetic 2024-T3 and Cu micro¬ 
electrode samples), arrays of Cu microelectrodes embedded in an 
insulator (CEI), and pure Cu. Unless otherwise noted, all electrodes 
were in the shape of a disk cut from sheets with diam 6.15 mm and 
area of 0.3 cm . The chemical composition (weight percent) of the 
2024-T3 material was as follows: 3.8-4.9 Cu, 1.2-1.8 Mg, 0.3-0.9 
Mn, 0.5 Si and Fe, 0.25 Zn, 0.1 Cr, and 0.05 Ti. The chemical 
composition (weight percent) of the solid solution alloys was as 
follows: Al-1% Cu-0.98 Cu, 0.001 Si, 0.002 Ti, 0.001 Zn; Al-4% 
Cu-3.96 Cu, 0.001 Si, 0.002 Ti, and 0.001 Zn. 

All sample disks were mechanically polished down to 3 |xm with 
successively finer Buehler Carbimet grinding paper that was lubri¬ 
cated with ethylene glycol. Further polishing was performed using a 
nonaqueous alumina slurry (suspension of Buehler Micropolish II 
deagglomerated alumina polishing powder in ethylene glycol) down 
to 0.05 pm. Following polishing, the sample disks were rinsed with 
ethanol. Water was not used during the polishing procedure owing to 
previous literature that reported accumulation of Cu on the surface 
and grooving around intermetallic particles during polishing of Al 
2024-T3 that used water in the slurry. 6 ' 9 

Synthetic 2024-T3 and Cu microelectrode samples. —After pol¬ 
ishing, some of the Al-1 wt % Cu disks were transferred to a clean- 
room where standard photolithographic techniques were used to 
make prescribed arrays of copper dots on the surface of the Al-Cu 
solid solution alloy disks and copper dots embedded in an insulator 


/ high pH zone Al 14(>l 

" fdi+4HX)44e~>40H’ 

dealloyed S-phase dealloyed matrix 



Figure 3. Schematic illustration of matrix dealloying in the vicinity of an 
S-phase particle on the surface of Al 2024-T3. The dealloyed S-phase serves 
as cathode active to oxygen reduction and the surrounding matrix undergoes 
Al dealloying causing a local enrichment of copper. 
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Figure 4. A schematic illustration of the growth zone of elevated pH around 
the dealloyed S-phase particle illustrating the planar diffusion front at times, 
r 0 , and the hemispherical shape of the zone at later times, / 3 . 5 is the 
thickness of the diffusion boundary layer. The shaded colors correspond to 
the shape and extent of the pH 9.5 zone as time evolves. 


(CEI). 14 Three different types of square lattice arrays of 100 nm 
thick copper dots were produced with dot diam of 1.8, 80, and 1000 
pun. Samples with 80 and 1000 pun diam Cu particles were made 
with square shaped matrix samples, 3 cm on a side. Since the dots 
covered an ~0.03 area fraction of the Al-1 wt % Cu alloy surface, 
the lattice sizes (interdot spacing) were 9.5, 400, and 5000 pun, 
respectively. 

Figure 6 (r = 0 h) shows a synthetic 2024 alloy made using an 
Al-1 wt % Cu solid-solution alloy with 1.8 pun diam Cu dots in a 
square array with nearest-neighbor spacing of 9.5 pun. 

RDE studies .—RDE studies were performed on the synthetic and 
real 2024 alloy samples in order to ascertain the amount of accumu¬ 
lated Cu on the surface as a function of rotation rate and time, and to 
measure oxygen diffusion-limited current densities as a function of 
rotation rate. Oxygen diffusion-limited current densities were also 
measured on planar Cu disks and CEI samples as a function of 
rotation rate. All of the experiments designed to measure oxygen 
diffusion-limited current densities were done in a 0.1 M Na 2 S0 4 
electrolyte under a prescribed rotation rate by performing an anodic 



t=0 h 



t=7h 


(b) 


Figure 5. Schematic illustration of the growth and overlap of zones of el¬ 
evated pH as a function of the cathode particle size, cl and spacing ,y. 5 is the 
thickness of the diffusion boundary layer. The shaded colors correspond to 
the shape and extent of the pH 9.5 zone as time evolves, (a) d ~ 8 , .v 
> 5, Particle separation is larger than the diffusion boundary layer thickness 
and zones of increasing pH are cut off by convection. Only the matrix im¬ 
mediately adjacent to the particles is exposed to a high enough pH for alka¬ 
line dissolution of Al. (b) d < 5, s < 8, Particle size and separation are 
smaller than the diffusion boundary layer thickness, the system behaves as a 
microelectrode array, causing the entire surface to activate to matrix dealloy¬ 
ing. 


t=20 h 

Figure 6. A time sequence of SEMs showing the morphology of a synthetic 
Al 2024 sample (1.8 pun Cu particles) after exposure to a 0.5 M NaCl 
solution. 
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sweep of the potential from the onset of hydrogen evolution at a 
sweep rate of 2 mV s” 1 . A saturated calomel electrode (SCE) was 
used as reference electrode and Pt wire as a counter electrode. A 
standard RDE device (Pine Instruments) together with a Pine Instru¬ 
ments bipotentiostat (model AFRDE5) was used. Corrosion treat¬ 
ments (A1 alloy 2024, Al-1 wt % Cu, synthetic 2024) at prescribed 
rotation rates and times were performed under free corrosion condi¬ 
tions in a 0.5 M NaCl solution. Afterward, the samples were taken 
out of the electrolyte, rinsed thoroughly with Barnstead Nanopure 
(>18 Mfl) water and transferred into another cell used for Cu assay 
using a lead underpotential deposition (Pb UPD) technique we pre¬ 
viously developed. 5 Typically, following Cu assay by Pb UPD, the 
samples were examined using both scanning electron and optical 
microscopy. 

RRDE studies. —RRDE studies were used to ascertain the opera¬ 
tion of the Cu dissolution/back-plating mechanism. The same Pine 
device together with the same reference and counter electrode was 
used as described in the section on RDE studies. A silver collection 
ring (area 0.125 cm 2 ) maintained at -0.80 V (SCE) was used to 
capture copper cations during free corrosion of samples (A1 2024- 
T3, Al-1 wt % Cu matrix, or synthetic 2024) in 0.5 M NaCl solution. 
Silver was chosen as the ring material instead of platinum in order to 
avoid the complicated stripping of Cu from Pt owing to the UPD 
phenomena, 16 as well as to minimize the hydrogen evolution. Each 
experiment involved immersion of the sample disk at a prescribed 
rotation rate and time. After the corrosion treatment, the ring along 
with the holder were rinsed with Nanopure water and transferred 
into a new cell containing 0.01 M HCI0 4 , where a potential of -0.4 
V (SCE) was immediately applied to the silver ring. In this cell the 
Ag ring served as the working electrode in a standard three- 
electrode configuration (SCE reference electrode, Pt counter elec¬ 
trode). After a steady-state current was attained, an anodic sweep 
was performed at a sweep rate of 20 mV s'" 1 . A positive anodic 
sweep limit of 0.2 V (SCE) was chosen in order to allow us to 
ascertain the onset of silver dissolution from the ring. The copper 
stripping wave typically appeared at potentials of 0.03-0.06 V 
(SCE). The advantage of introducing a second electrochemical cell 
was to avoid the formation of AgCl at potentials close to that of 
copper stripping. 

All the solutions were made with high purity grade chemicals 
and Nanopure water. A BAS CV-27 model potentiostat was used for 
the analytical electrochemical experiments (Cu stripping from the 
Ag ring and Cu assay by Pb UPD) and a Nicolet 310 digital oscil¬ 
loscope was used to store data. All glassware used for experimenta¬ 
tion and the preparation of electrolytes was cleaned with concen¬ 
trated HN0 3 and H 2 S0 4 at 70°C, then rinsed using doubly distilled 
and Nanopure water. 


Results and Discussion 

Figure 7 shows results for the Cu content determined by Pb UPD 
assay on synthetic and real 2024 samples as a function of treatment 
time in 0.5 M NaCl at zero rotation rate. We note that the assay was 
accurate to within ±0.02 area fraction of surface coverage. 15 These 
results pose the major question addressed in this work, namely, what 
are the causes of the differences in the time-dependent surface cov¬ 
erage of Cu between the synthetic and the real alloy? The following 
is a list of the possible causes of this discrepancy. 

Dissolution and back plating contribution: porous S-phase. —As 
discussed in the first part of this paper, a requirement of the disso¬ 
lution and back-plating of copper redistribution is the mechanical/ 
electrical detachment of copper particles. This allows the detached 
copper particle to adopt the corrosion potential of copper in the 
electrolyte leading to dissolution of copper and back-plating of Cu 2+ 
ions from solution. There are many anodic processes available to 
support back-plating, including second-phase and/or matrix dealloy¬ 
ing, matrix oxidation, and double-layer processes. One likely source 
of copper liberation from A1 2024-T3 that is unavailable to the syn¬ 
thetic alloy is dealloyed S-phase. In general, dealloying will not 



Figure 7. Cu surface coverage, determined by Pb UPD assay, on synthetic 
and real A1 2024 samples as a function of exposure time in 0.5 M NaCl. 


always produce a mechanically stable bicontinuous structure. For 
example, consider the dealloying of a Ag {)9 Au (U alloy. Since the Au 
atom fraction is below the 3-D site percolation threshold (corre¬ 
sponding to 20 atom % gold), the formation of a mechanically stable 
bicontinuous gold structure at this alloy composition is not possible. 
Instead, dealloying of silver produces 3-D gold clusters that become 
disconnected from the corroding substrate. The Al 2 CuMg, S-phase, 
contains 25 atom % copper, and has an orthorhombic crystal struc¬ 
ture. From a percolation point of view, this intermetallic is a border¬ 
line case in the sense of resulting in mechanically stable bicontinu¬ 
ous dealloyed structures. During S-phase dealloying, we cannot rule 
out the possibility that nanometer-scale ligaments of copper become 
disconnected from the Cu-rich porous sponge evolving during the 
dealloying process. As discussed below, another possible source of 
liberated copper for both the real and synthetic 2024 is matrix deal¬ 
loying. In order to evaluate the relative contribution of these pro¬ 
cesses to copper redistribution, we performed RDE and RRDE stud¬ 
ies. 


Surface activation. —Figure 7 shows that as the particle spacing 
in the synthetic alloy gets smaller the amount of copper on the 
surface increases. As illustrated in Fig. 5, there are two length scales 
that determine the behavior of the synthetic alloys: one is the sepa¬ 
ration between the particles on the surface (at fixed area fraction this 
sets the particle diameter) and the other is the thickness of the dif¬ 
fusion boundary layer. The requirement for the entire surface to 
“activate” to oxygen reduction is that the Cu particles on the surface 
must behave collectively as a microelectrode array. This behavior 
occurs when the Cu dot separation in the synthetic alloy is much less 
than the diffusion boundary layer thickness. In this case, the oxygen 
diffusion-limited current density of the Cu microelectrode array is 
equal to that of a planar copper electrode. If the interparticle sepa¬ 
ration is too large, then the oxygen reduction current density is low¬ 
ered, yielding corresponding reductions in the corrosion rate and 
subsequent copper accumulation on the surface of the synthetic al¬ 
loy. We have examined these effects by RDE measurements of oxy¬ 
gen diffusion-limited current densities on all five of the samples 
described in the Experimental section. An obvious route to test the 
effects described above is to simply produce model alloys with 
smaller Cu dot separations, while maintaining a 0.03 area fraction of 
Cu on the surface. Unfortunately, the synthetic alloy with 1.8 (xm Cu 
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Figure 8. Schematic illustration of the local pH probe experiment. 


dot diam already represents the order of magnitude limit of feature 
sizes and dimensions available via conventional lithography routes. 

Copper in solid solution .—The synthetic alloy could contain less 
Cu in solid solution then the real A1 2024. In principle, the amount 
of Cu in solid solution can be examined via high resolution analyti¬ 
cal transmission electron microscopy (TEM), however in practice 
these measurements are difficult. This hypothesis was examined in 
ancillary experiments, using synthetic alloys containing 4 wt % Cu 
in solid solution. These results are presented at the end of this paper. 

Ancillary issues. —Prior to presenting the main body of our ex¬ 
perimental results, it was necessary to consider three additional is¬ 
sues. The first issue concerned the behavior of the synthetic alloy 
vis-a-vis retention of the copper dots on the surface, during immer¬ 
sion in the chloride electrolyte. 

Figure 6 shows a time sequence of scanning electron micro¬ 
graphs (SEMs) of the surface of the synthetic alloy following im¬ 
mersion under free corrosion in the chloride electrolyte. For immer¬ 
sion times of order a few hours, there is little damage to the surface 
and the array of Cu particles stays intact. While periphery grooving 
around the copper dots is apparent in samples that were immersed 
for times of order 7 h or more, close inspection reveals that each of 
the “pits” retains a copper dot at the bottom. We would expect the 
copper dots to be undercut eventually by corrosion. Figure 7 shows 
that the relevant time scale to achieve maximum copper surface 
coverage in a stagnant electrolyte is of the order of 100 min, and 
over these time scales, the synthetic alloys retain integrity. 

Another issue relates to the steady-state pH that evolves within 
the diffusion boundary layer during corrosion of A1 2024-T3 in a 
nominally neutral electrolyte. In order for the matrix dealloying 
mechanism to operate, the pH in the diffusion boundary layer should 
be greater than —8.0 for aluminum to dissolve via the soluble 
A10 2 ” species. 17 We measured the local pH in the vicinity of 
S-phase particles on the surface of A1 alloy 2024-T3 by using the tip 
of a scanning tunneling microscope (STM) as an in situ local pH 
probe as illustrated in Fig. 8. A 99.5% pure Johnson-Matthey Pt- 
20% Ir wire (0.25 mm diam) was used as a working electrode and a 
Pt wire (1 mm diam) served as the counter electrode. In all the 
experiments, the potential was measured vs. SCE. Calibration solu¬ 
tions were made using 0.1 M NaCl + x M NaOH where x was 
chosen to provide us with pH 6.5, 7.5, 8.5, 9.0, 10.0, 10.5, 11.2, and 
12.8. The pH of 0.1 M NaCl (x = 0) was 5.3. In these experiments, 
the potential was swept from 300 to 1500 mV (SCE) at a sweep rate 
of 20 mV s' 1 . 

Figure 9 represents the current-voltage (I-V) behavior of a Pt- 
20% Ir wire in 0.1 M NaCl + x M NaOH. The first curve where no 
NaOH was added (pH 5.3) demonstrates a flat double layer region 
with the onset and further development of an oxygen evolution pro¬ 
cess. All other curves reveal a current peak at potentials negative to 
the potential of the oxygen evolution process. Similar peaks ob¬ 
tained on polycrystalline and monocrystalline platinum electrodes in 
alkaline solution by others 18 have been related to oxygen adsorption 
that precedes the oxygen evolution. The peak current densities in 
Fig. 9 are plotted vs. pH in Fig. 10. A simple linear regression 



Figure 9. I-V behavior of a Pt-20% Ir microelectrode in 0.5 M NaCl 
+ x M NaOH. Sweep rate, 20 mV s" 1 . 


analysis was used to fit the data that served as a standard for the in 
situ STM measurements. 

In situ STM measurements over and immediately adjacent to 
S-phase particles in the surface of A1 alloy 2024-T3 immersed in a 
neutral chloride electrolyte revealed that the pH tended toward a 
steady-state value of —9.5 within 20 min of immersion. 

The third issue is concerned with the role of second phase par¬ 
ticles ( e.g ., FeMnCuAI, Al 2 Cu, etc.) other than the S-phase in serv¬ 
ing as a source of redistributed copper in real 2024-T3. In ancillary 
studies, Ford et a/. 19 used analytical TEM to characterize the 
2024-T3 material used in this study following corrosion treatments 
in a neutral chloride solution in order to determine the reactivity 
(, i.e ., the ability of the particle to serve as a source of redistributed 
Cu) of the S-phase and other constituent particles such as those 



pH 


Figure 10. Peak current density of a Pt-20% Ir microelectrode in 0.5 M 
NaCl + x M MNaOH as a function of pH. 
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Figure 11. Cu surface coverage (Pb UPD assay) on real A1 2024-T3 after 
free corrosion in 0.5 M NaCl as a function of exposure time at the indicated 
rotation rates. 


containing FeMnCuAl. They found that the S-phase dealloyed and 
that FeMnCuAl dispersoids were stable. Phases such as Al 2 Cu, even 
if dealloyed will produce a copper-rich remnant that would remain 
mechanically intact since the copper content is well above the per¬ 
colation threshold (see the Results and Discussion section on disso¬ 
lution and back-plating contribution: porous S-phase). Mazurk- 
iewicz and Piotrowski 20 examined the behavior of monolithic Al 2 Cu 
in 0.5 M S 04 - electrolyte in an RRDE configuration. By observing 
color change, they reported large quantities of Cu captured on the 
ring indicating the operation of a Cu back-plating mechanism. The 
implication is that the Al 2 Cu phase may be an important source of 
redistributed copper in Al-Cu alloys. Our own RDE and RRDE re¬ 
sults presented next in the sections on RDE studies, and RRDE 
studies: copper redistribution contradict this conclusion. 11 Therefore, 
we consider the S-phase as the only particle constituent phase that 
could serve as a significant source of redistributed copper. 

RDE studies. —Figure 11 shows results for copper accumulation 
on the surface of A1 2024-T3 as a function of free corrosion immer¬ 
sion time (0.5 M NaCl) and rotation rate. The data joined by solid 
lines were obtained by using the same sample for a prescribed time 
and rotation rate. Following copper assay by Pb UPD, the sample 
was washed using procedures described in the Experimental section 
and was directly reimmersed at the same rotation rate. The time axis 
represents cumulative exposure times for this data set. The data set 
joined by dashed lines were obtained from samples that were im¬ 
mersed for a prescribed time and rotation rate, subjected to copper 
assay by Pb UPD and repolished using procedures described in the 
Experimental section. The time axis represents total exposure times 
for this data set. The results for copper accumulation are insensitive 
to the details of the methods used to obtain the data sets. The results 
indicate that for immersion times less than about 10 min the accu¬ 
mulated copper increases with rotation rate. For immersion times 
greater than about 10 min the accumulated copper increases with 
rotation rate in the range of zero to 400 rpm and then either saturates 
or shows a slight decrease. The coverage saturates at an equivalent 
monolayer of copper on the surface. 

d We suspect that their results were affected significantly by the manner in which the 
A1 2 Cu disk was embedded into an insulator resulting in large compressive radial 
stresses that caused the remnant porous Cu to break apart during dealloying of Al. 


The general effect of rotation rate on copper coverage is easily 
understood. Since the thickness of the diffusion boundary layer de¬ 
creases as the reciprocal of the square root of rotation rate, the 
oxygen diffusion-limited current density increases. Thus, the corro¬ 
sion rate of the sample increases with rotation rate and we can 
expect the surface copper coverage to increase accordingly. 

Figure 12 shows optical micrographs of the Al 2024-T3 surface 
following 60 min of free corrosion in the 0.5 M NaCl electrolyte at 
the indicated rotation rates. In the stagnant electrolyte, the vast ma¬ 
jority of corrosion damage is localized in the vicinity of S-phase 
intermetallic particles imbedded in the alloy matrix. Since the cor¬ 
rosion rate increases with increasing rotation rate, the corrosion 
damage should increase accordingly. We observe this behavior up to 
rotation rates of about 400 rpm where qualitatively the corrosion 
damage is a maximum as accessed by the overall degree attack on 
the alloy matrix between the particles. At higher rotation rates, the 
corrosion damage begins to localize around the particles and at 5000 
rpm the damage is mainly confined to the region immediately sur¬ 
rounding the particle. This behavior is a result of two competing 
effects associated with the increased rotation rate or decreasing dif¬ 
fusion boundary layer thickness. As the diffusion boundary layer 
decreases, the oxygen diffusion-limited current density increases, as 
does the corrosion rate. Eventually, as the rotation rate increases, the 
diffusion boundary layer thickness gets to be of order the size of the 
cathode (particle) separation on the surface of the alloy. The hemi¬ 
spherical regions of oxygen diffusion (which mirror the zones of 
increasing pH) are thus terminated. As a result, the cathodic current 
density is reduced below the mass-transport limited rate. Concur¬ 
rently, the hemispheres of increasing pH, emanating from the par¬ 
ticles (see Fig. 4), are cut off by convective flow and the matrix 
regions separating the particles never achieve the high pH required 
for matrix dissolution. Owing to these competing effects, beyond a 
certain rotation rate, the area fraction of the surface that shows sig¬ 
nificant dissolution or corrosion damage decreases with increasing 
rotation rate. Thus we expect to see the maximum corrosion damage 
in the matrix, between the particles active to oxygen reduction, at 
some intermediate rotation rate. The rotation rate corresponds to a 
diffusion boundary layer thickness defining an average separation 
between active particles in Al 2024-T3. Application of the Levich 
equation indicates that this length scale is of order 10-20 pan. In 
comparing conclusions based on the results presented in Fig. 11 and 
12, there is an unresolved issue. Namely, why is it that the copper 
accumulation increases with increasing rotation rate, while the sur¬ 
face corrosion damage saturates in the vicinity of 400 revolutions 
per minute (rpm)? We return to this question in a later section of this 
paper. 

Figure 13 shows results for copper accumulation on the surface 
of the synthetic 2024 alloy (9.5 p,m Cu dot spacing) as a function of 
free corrosion immersion time and rotation rate. The time- 
dependent, accumulated copper is independent of rotation speed up 
to 200 rpm and then decreases with increasing rotation speed. This 
behavior is a direct consequence of our previous discussion regard¬ 
ing Fig. 12. At 1500 rpm, the diffusion boundary layer thickness is 
about 13 |xm, which is close enough to the interdot distance so that 
the zones of increasing pH are cut off by convection. Under identi¬ 
cal conditions, the area fraction of accumulated copper on the sur¬ 
face of the synthetic 2024 alloy is consistently lower that that ob¬ 
served on the real 2024. We note that if phases such as AI 2 Cu 
(which likely exist in both the synthetic and real 2024) were impor¬ 
tant in the copper redistribution process, then the results presented in 
Fig. 11 and 13 would have been quantitatively similar. Figure 14 
shows optical micrographs of the synthetic 2024 following 60 min 
of free corrosion in the 0.5 M NaCl electrolyte at the indicated 
rotation rates. In general, one observes that the corrosion becomes 
more localized with increasing rotation rate. 

RRDE studies: copper redistribution. —The results presented 
above demonstrated an influence of the rotation rate on the copper 
redistribution process. RRDE experiments were performed in order 
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Figure 12. Optical micrographs of real AI 2024 samples after treatment in 0.5 M NaCl for 60 min at the indicated rotation rates. 


to determine if the Cu dissolution and back-plating mechanism is 
operative in the Cu redistribution process. In order to ascertain this, 
a fixed potential was maintained on the silver ring that was more 
negative than the corrosion potential of the Al matrix. The applied 
potential allowed the silver ring to enter into competition with the Al 
matrix in the Cu back-plating process. A relatively long corrosion 



Figure 13. Cu surface coverage (Pb UPD assay) on synthetic Al 2024 
samples following free corrosion in 0.5 M NaCl as a function of exposure 
time at the indicated rotation rates. 


time (90 min) was chosen in order to insure that if there were any 
Cu 2+ ions in solution a detectable amount would be deposited on the 
Ag ring. A back of the envelope estimate shows that at least 50 
equivalent monolayers of Al should be dissolved from the synthetic 
2024-T3 in order to register one Cu monolayer on the ring surface. 

Figure 15 shows results of Ag-ring anodic stripping waves ob¬ 
tained in 0.01 M HC10 4 following open circuit exposure of a pure 
copper disk and Al alloy 2024-T3 disk in 0.5 M NaCl at 1000 rpm. 
The I-V behavior of a bulk polycrystalline copper electrode is also 
shown. A comparison of these curves allowed us to conclude that the 
stripping wave associated with corrosion of the Al 2024-T3 resulted 
from copper dissolution. The charge associated with the stripping 
wave of the Al 2024-T3 was about 2% of that obtained from the 
pure Cu disk. Figure 16 is a set of copper stripping curves from the 
Ag ring (obtained in 0.01 M HC10 4 electrolyte) after free corrosion 
of Al 2024-T3 and Al-1 wt % Cu (matrix) disks in 0.5 M NaCl for 
90 min. The amount of copper captured on the silver ring during 
corrosion of both the Al 2024-T3 and Al-1 wt % Cu matrix samples 
in a stagnant 0.5 M NaCl electrolyte was within the detection limits 
of the technique. Similar results were obtained by Buchheit et al. 2{ 
using a Pt ring. The amount of Cu captured by the Ag ring was 
found to be a strong function of rotation rate. Figure 17 obtained by 
integration of the stripping waves of Fig. 16 shows an almost linear 
increase of the captured copper in the range of 1000-2500 rpm re¬ 
gion. Identical results to those obtained for the Al-1 wt % Cu matrix 
disks were also obtained for the synthetic alloys, reaffirming that the 
only source of Cu available for redistribution from the synthetic 
alloy is copper in solid solution. These results demonstrate that the 
S-phase on the surface of Al 2024-T3 is an important source of 
redistributed copper via the dissolution and back-plating mecha¬ 
nism. The Al-1 wt % Cu matrix disks emitted relatively modest 
amounts of copper into the chloride electrolyte. Figures 18 and 19 
show the Cu stripping waves as a function of time at 1500 rpm and 
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Figure 14. Optical micrographs of synthetic A1 2024 samples after free corrosion in 0.5 M NaCl for 60 min at indicated rotation rates. 


Fig. 20 summarizes these results after integration of the stripping 
waves. Note, that for the real 2024, copper is present on the ring 
almost immediately, whereas it takes about 60 min for a detectable 
quantity to appear for the Al-1 wt % Cu-matrix samples or the 
synthetic alloy. The RRDE results (Fig. 16-20) are consistent with 
the RDE results in that if phases such as Al 2 Cu were important to 



Potential. mV fSCE] 

Figure 15. Copper stripping waves obtained from the Ag ring in a 0.01 M 
HCIO 4 solution. Copper on the Ag ring accumulated from free corrosion of 
a pure copper disk and a real A1 2024 disk at 1000 rpm in 0.5 M NaCl. The 
I-V behavior of a bulk copper electrode in a 0.01 M HCIO 4 solution is shown 
for comparison. 


the copper redistribution process, the results for the matrix alloy, the 
synthetic 2024, and the real 2024 would have been quantitatively 
similar. 

RDE studies: oxygen diffusion limited current 

densities. —Regardless of the operative copper distribution mecha¬ 
nisms, intermetallic particles serve as the primary sites for oxygen 
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Figure 16. Copper stripping curves from a silver ring after free corrosion of 
real 2024 and matrix samples for 90 min in 0.5 M NaCl at the indicated 
rotation rates. 
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0 50 100 150 200 250 300 Figure 19. Copper stripping curves from a silver ring after free corrosion of 

j -1 matrix samples in 0.5 M NaCl for different times at 1500 rpm. 


Figure 17. Current densities of copper stripped from the Ag ring in 0.01 M 
HC10 4 solution after 90 min of free corrosion in 0.5 M NaCl of matrix, 
synthetic, and real A1 2024 samples as a function of rotation rate. 


reduction on a corroding A1 2024-T3 surface. In this regard, we 
assume that dealloyed S-phase particles are most important. How¬ 
ever, our major conclusions are unaffected by contributions from 
other types of particles. 22,23 The oxygen reduction rate balances the 
rate at which any anodic process related to copper redistribution 
operates. Therefore, this parameter must be well characterized if one 
is to develop a quantitative model of corrosion and copper redistri¬ 
bution in Al-Cu-X alloys. Furthermore, RDE measurements of oxy¬ 
gen diffusion-limited current densities allow us to determine if the 
synthetic alloy geometry (1.8 jim diam Cu disks in a square array 
with 9.5 pim near-neighbor spacing) is behaving as an array of mi- 



Potential, mV (SCE) 


Figure 18. Copper stripping curves from a silver ring after free corrosion of 
real A1 2024 samples in 0.5 M NaCl for different times at 1500 rpm. 


croelectrodes. For example, at low enough rotation rates, the oxygen 
diffusion limited current density of CEI samples should be identical 
to that of a planar Cu electrode. We note that since both anodic and 
cathodic processes are operative in the synthetic 2024, the net ca¬ 
thodic current density will be less than that of the CEI or planar Cu 
samples. 

In a conventional four-electron reduction mechanism, under aer¬ 
ated conditions, the oxygen diffusion-limited current density in a 
stagnant electrolyte is a poorly defined quantity in the range of 20- 
100 \iA cm -2 . Additionally, recent studies in our laboratory demon¬ 
strate that the magnitude of the oxygen diffusion-limited current 
density on planar copper electrodes is not a function of electrolyte 
pH. 24 We found that measurements in pH 1.5-14 electrolytes yielded 
values between 20 and 30 jxA cm -2 . In these electrolytes, at a po- 



Time, min 

Figure 20. Charge densities as a function of time of stripped copper from 
the silver ring in 0.01 M HC10 4 solution after free corrosion of matrix and 
real A1 2024 samples in 0.5 M NaCl solution at 1500 rpm. 
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Figure 23. Cu surface coverage (Pb UPD assay) on real, synthetic A1 2024 
(made using a 4 wt % Cu matrix), and a 4 wt % Cu matrix sample as a 
function of exposure time in 0.5 M NaCl. 


The CEI sample shows distinctly different behavior. In general, 
for a microelectrode array, the limiting current will follow, i L 
~ L/X^) J * [*L.pianar electrode (^)J» where /(to) is a rotation rate- 
dependent shape function defining the volume fraction per unit time 
of the diffusion boundary layer supplying oxygen for the reduction 
process. The shape function equals 1.0 in the limit of low rotation 
rates and then decreases with increasing rotation rate. Shape func¬ 
tion effects are apparent in the optical micrographs (Fig. 12 and 14) 
of the corroded real and synthetic alloys. 

Conclusions 

The first issue (of the questions and issues posed by the results 
presented in Fig. 6) that we address is that of “surface activation” of 
the synthetic alloy (Al-1 wt % Cu) composed of 1.8 p,m diam Cu 
disks in a square array with 9.5 p,m near-neighbor spacing. Our 
studies of the oxygen diffusion-limited current density demonstrate 
that this geometry behaves as a microelectrode array in the limit of 
a stagnant or very weakly stirred electrolyte. Under this constraint, 
the microelectrode array provides the same oxygen diffusion-limited 
current density as a planar Cu electrode. Thus, in a stagnant electro¬ 
lyte, this synthetic 2024 alloy is as efficient as the real 2024 for 
oxygen reduction. 

Our RRDE results and those of Buchheit et al , 21 demonstrate 
the operation of a dissolution/back-plating copper redistribution 
mechanism during corrosion of Al 2024-T3. This process is trig¬ 
gered by very small copper particles of order tens of nanometers in 
size that become mechanically disconnected from dealloyed S-phase 
intermetallic particles. As long as the remnant dealloyed copper 
sponge stays mechanically and electrically connected to the sub¬ 
strate, no copper dissolution can take place. In this case, we can 
disregard particle size effects on dissolution since the corrosion po¬ 
tential is almost 1 V less than the reversible potential of copper. 10,27 

The RDE studies of the oxygen diffusion-limited current density 
demonstrate that the synthetic 2024 (1.8 pm diam Cu disks in a 
square array with 9.5 pm near-neighbor spacing) geometry behaves 
as a microelectrode array for oxygen reduction. Furthermore Fig. 22 
shows that independent of rotation rate, oxygen reduction on the 
.surface of the real and synthetic alloys is quantitatively similar. Dur¬ 
ing immersion of these alloys in 0.5 M NaCl, the only physical 
process that we can identify taking place on the surface of the real 


2024 that does not occur on the synthetic 2024 is the operation of 
the dissolution/back-plating mechanism of copper redistribution. 

We address the final point relating to the amount of copper in 
solid solution in the real and synthetic 2024. Since Al 2024-T3 is 
naturally aged, no single well-defined number can characterize the 
level of copper in solid solution. A large fraction of the ~4 wt % Cu 
in this alloy eventually is taken up by coherent or incoherent second 
phases. As noted earlier, the literature indicates that a reasonable 
value for the copper in solid solution is 0.2-0.5 wt %? Thus, we 
believe that our synthetic alloy is a reasonable model in this regard. 
However, it is easily demonstrated as shown in Fig. 23, that as the 
copper content in solid solution of a synthetic alloy increases to 4 wt 
% the redistributed copper on the surface of the synthetic alloy 
closely approximates the behavior of the real 2024. In this case, the 
matrix dealloying mechanism alone is sufficient to quantitatively 
account for the copper redistribution process. On-going work in our 
laboratory using quantitative analytical high resolution TEM is 
aimed at ascertaining the amount of copper in solid solution of the 
Al 2024-T3 alloys that we have been using. To date, we find that the 
amount of Cu in solid solution is spatially heterogeneous and ranges 
from just below 1 to 4 wt % and so presently we cannot make any 
definitive conclusions on this issue. 28 

Finally, the differences in behavior for the surface coverage of 
Cu between 2024-T3 and the synthetic 2024 alloy is attributed to the 
operation of a dissolution/back-plating mechanism triggered by me¬ 
chanical disconnection of Cu particles from dealloyed porous 
S-phase. Comparisons of the behavior of these two alloys allow us 
to conclude that the matrix dealloying and dissolution/back-plating 
contribute about equally to the copper redistribution process during 
corrosion of aluminum alloy 2024-T3 in 0.5 M NaCl. 
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Abstract 

There are many important processes in corrosion for which the diffusion -limited 
current density of oxygen reduction, i L , plays a dominant role in terms of kinetic control. 
The conventionally believed magnitude of i L (for the four-electron reduction mechanism) 
which can be found in many corrosion textbooks is in the range of 50-100 |iA cm ', but 
the origins of this range of values are a bit mysterious. Previous research in our group 
aimed at ascertaining i L (under stagnant conditions) on a planar Cu electrode and a Cu 
microelectrode array in a naturally aerated 0.1M Na 2 S 0 4 electrolyte found that it was in 
the range of 20 - 30 |iA cm' 2 . Obvious reasons for this discrepancy may involve a two - 
electron mechanism in the reduction or the presence of a surface film (an oxide) that 
retards the kinetics of the diffusion - limited reduction process. The work reported on 
herein was aimed at sorting out these issues. In situ STM was used to characterize the Cu 
surface at relevant potentials. RDE studies were used to measure i L in a naturally aerated 
0.1M Na 2 S0 4 electrolyte as a function pH for both Pt and Cu electrodes. By comparing 
results for Pt and Cu we conclude that oxygen reduction occurs on a Cu surface via the 
four - electron mechanism. We found that while there is an underpotentially deposited 
oxygen adlayer on the Cu surface for pH values of ~ 9.5, this had no effect on i L . The 



oxygen diffusion - limited current density was found to be independent of pH (in the 
range 1.5 - 14). Finally, we conclude that our previous determination of k in a stagnant 
electrolyte actually reflect an accurate range of values of the oxygen diffusion - limited 
current density. 
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Introduction 


Oxygen reduction is important in corrosion since the kinetics of this process often 
determines general corrosion rates and damage evolution (e.g., pitting, crevice corrosion, 
stress -corrosion, corrosion fatigue) in a number of commercially important alloys. 
Nevertheless, as briefly discussed below, there have been relatively few well - conceived 
experiments aimed at measuring the rate of oxygen reduction under conditions directly 
relevant to corrosion phenomena. To date, most of the significant work in this arena was 
aimed at studying the catalytic activity of Pt (hkl) to the oxygen reduction process using a 
rotating ring - disk electrode (RRDE) configuration under oxygen saturated conditions 
[ 1 ]. 


A significant motivation for the study that we report on herein is the key role of 
oxygen reduction in the corrosion of Aluminum Alloy 2024-T3. In this alloy, the 
corrosion rate is controlled by oxygen reduction on spatially separated microscopic 
intermetallic particulate phases that are highly enriched in copper owing to de-alloying 
[2-7]. On immersion in a corrosive environment these intermetallic particles (e.g., S- 
phase A^CuMg [2,4]) support oxygen reduction, presumably according to the reaction 

2 H 2 O + O 2 + 4e~ = 4(OH), which results in a pH increase in the neighborhood of the 
local micro - cathodes. As the pH increases to about 9-9.5 [8-10], the passive oxide on 
the surface of the alloy matrix chemically dissolves causing the aluminum matrix to 

dissolve via the soluble A10 2 anion. We note that the corrosion potential of A1 alloy 
2024-T3 in a 0.5 M chloride or sulfate electrolyte is -800 mV to -600 mV (SCE) [11]. 


3 



The general corrosion rate of the alloy and the subsequent extent of damage will depend 
on kinetics of the oxygen reduction at the copper enriched particles. In a recent study [2], 
we compared the kinetics of oxygen reduction on A1 2024 T3, a synthetic model A1 2024 
T3 alloy, a planar Cu electrode and a lattice of copper particles embedded in an insulator 
(CEI samples) in 0.1 M Na 2 S 04 . Both the synthetic alloy and the CEI sample contained 
1.8 pm diameter Cu disks arranged in a square lattice with 9.5 pm near-neighbor spacing. 
This particle size and density was chosen to simulate the surface coverage of the S-phase 
on A1 Alloy 2024-T3 [2], We found that this configuration of Cu particles behaved as a 
microelectrode array for oxygen reduction. The oxygen diffusion - limited current 
density */, on the CEI microelectrode array and on a pure planar copper electrode in a 
stagnant, naturally aerated electrolyte was 20-30 pA cm' 2 [2]. These numbers differ 
significantly (by a factor of 4-5) from the numbers commonly accepted and used in the 
corrosion literature (i.e., 50 - 100 pA cm' 2 [12]) and the immediate question is why? One 
possibility is that in a pH ~ 9 electrolyte in the near electrode vicinity, resulting from the 
ongoing oxygen reduction, a Cu surface may be oxidized so oxygen reduction may be 
under combined mass-transport and surface reaction - limited control. 

In a truly stagnant electrolyte the diffusion boundary layer thickness (8) tends to 
infinity and it would asymptotically approach a negligible value. However, in a real 
system ii is never zero because of naturally occurring convection and a finite 8 (albeit 
poorly defined) results. 
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Existing results for oxygen reduction on copper were obtained in oxygen saturated 
solutions [13,14]. They suggest a four-electron pathway for the process as well as an 
influence of the surface oxidation state on the reaction kinetics [13]. An influence on the 
reaction kinetics on pH as well as anion was also reported [14]. 

The study reported on herein is aimed at investigating the oxygen reduction reaction, 
taking place on a pure copper surface in naturally aerated 0.1 M NaiSClj solution as a 
function of pH. Separate experiments performed on platinum under identical 
experimental conditions allowed for a comparison of the two metals as substrates for the 
oxygen reduction reaction. The Rotation Disk Electrode (RDE) technique was used in 
this investigation. Additionally, electrochemical scanning tunneling microscopy 
(ECSTM) experiments were used to characterize the morphological changes on a Cu 
(111) surface at potentials corresponding to the corrosion potential of A1 2024-T3. 

2. Experimental Methods 

2.1 Samples 

A Cu (111) single crystal 2 mm thick and 10 mm in diameter was used in our 
ECSTM study. The crystal was electrochemically polished in 5:4:1 H 3 PO 4 : water : ethyl 
alcohol solution prior to each ECSTM experiment. Anodic DC current (0.1 - 0.15 A cm' 
2 ) was applied for 20 - 30 seconds. After the electrochemical polishing, the crystal was 
thoroughly rinsed with Barnstead Nanopure® (>18 MC2) water and covered with a 
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droplet of pH 2 sulfuric acid in order to protect it from oxidation. Polycrystalline Cu and 
Pt disks with diameters 0.615 cm and areas of 0.3 cm 2 were cut from sheets and served as 
working electrodes in the RDE experiments. The disks were mechanically polished down 
to 3 pm on BUEHLER Carbimet® grinding paper that was lubricated with ethylene 
glycol. An alumina slurry (suspension of BUEHLER MICROPOLISH II 
Deagglomerated alumina polishing powder 0.05 pm in ethylene glycol), was used as the 
final mechanical polish. Following polishing the sample disks were rinsed with ethanol. 
The platinum disk was hydrogen flame-annealed for 5 minutes. 

2.2 Rotating Disk Electrode Studies 

RDE studies were performed on polycrystalline Cu samples in order to measure 
oxygen diffusion-limited current densities as a function of rotation rate and pH. The pH 
employed was 1.5, 3, 5.5, 7.5, 10, 12, and 14. Solutions with pH values of 1.5 and 3 were 
made using 0.1 M Na 2 S0 4 acidified with H 2 S0 4 while solutions with pH values of 7.5; 

10; 12 and 14 were made using 0.1 M Na 2 S0 4 adjusted with NaOH. RDE experiments 
were performed also on a Pt polycrystalline electrode in 0.1 M Na 2 S0 4 (pH 5.5) and in 
0.1 M Na 2 S0 4 adjusted to pH 10. Each of the experiments designed to measure oxygen 
diffusion-limited current densities were at constant pH under a prescribed rotation rate by 
performing separately a cathodic as well as an anodic sweep at a sweep rate of 2 mV s' 1 . 
The cathodic sweeps were performed starting from open circuit potential while the anodic 
sweeps had an initial potential just above the onset of hydrogen evolution. A saturated 
calomel electrode (SCE) was used as a reference electrode and a Pt sheet as a counter 
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electrode. All the potentials in the text are reported vs SCE unless stated otherwise. A 
standard RDE device (Pine Instruments) together with a Pine Instruments bi-potentiostat 
(model AFRDE5) were used and a Nicolet 310 Digital Oscilloscope served to store data. 

2.4 ECSTM Experiments. 

The ECSTM experiments were carried out using a Molecular Imaging Pico STM with 
a 300S scanner and a Molecular Imaging model 300S Pico Bipotentiostat. The cell was 
made out of Teflon and had a volume of ca. 1.5 cm' exposed an area of ca. 0.3 cm . The 
tunnelling tips were made of a polycrystalline Pt-20%Ir wire electrochemically etched in 
a solution of 1:3 saturated CaCl 2 : water at about 21 V (AC). The tips were coated with 
Apiezon wax, leaving only the very end of the tip exposed to the electrolyte, in order to 
reduce the background current to numbers lower then 20 pA. 

The working electrode in the ECSTM study was a Cu (111) single crystal. Platinum 
wires served as a reference and counter electrodes. They were cleaned with hot sulphuric 
and nitric acids and hydrogen flame annealed prior to the experiment. The reference 
electrode had a relatively stable potential of 50±20 mV. The electrolyte used was 0.1 M 
Na 2 S 04 adjusted to pH 10. In the ECSTM experiments, a potential of-1250 mV was 
applied immediately after introducing the solution to the cell, thus preventing the surface 
from spontaneous oxidation. There were two kinds of experiments. The first one was an 
examination of the Cu surface by increasing the potential in step-wise increments from - 
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1100 mV to -750 mV. In the second one we kept the potential at -800 mV and followed 
the morphology evolution with time. 

All the solutions were made with high purity grade chemicals and Barnstead 
Nanopure® water. All glassware and Teflon ECSTM electrochemical cell used for 
experimentation and the preparation of electrolytes were cleaned with preheated (70°C) 
concentrated HNO 3 and H 2 SO 4 and then rinsed using doubly distilled and Barnstead 
Nanopure® water. 

3. Results and Discussion 

STM Results 

Figure 1 shows a set of STM images of a Cu (111) surface in 0.1 M Na 2 S 04 adjusted 
to pH 10 solution recorded while increasing the potential stepwise from -1100 mV to - 
750 mV. The first two images (la and lb) represent a bare copper surface that consists 
mostly of monoatomic steps separating atomically flat terraces. No changes could be 
detected on the surface while increasing the potential from -1100 mV to -1000 mV. The 
first hint of a change was detected at a potential of-900 mV (image lc). There was a 
slight increase in size of the topmost cluster in the middle portion of the image. More 
pronounced differences were found on image Id where at a potential of-800 mV all 
terraces have grown out in a lateral direction. The formation of an underpotentially 
deposited (UPD) oxide layer was observed at a potential of-750 mV (image le) where 
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the entire surface was covered by high density oxide clusters the height of which 
correlates with the size of the Cu-0 dimmer. A picture of the reduced Cu (111) surface 
obtained by switching the potential to -1300 mV is included at the end of this set of 
images. Similar results were obtained by Maurice et al. [15] on Cu (111) surfaces in 0.1 
M NaOH at potentials of-950 mV to -900 mV. In their work the authors explain the 
lateral growth of the terraces as a result of “accumulation, at the step edges, of Cu atoms 
ejected from the adlayer superstructure due to reconstruction of the outermost Cu plane 
of the substrate”. It was also suggested that “the same mechanism is responsible at the 
end of the adsorption process for the formation of Cu adislands on terraces when the 
adlayer grows in adsorbate-free areas surrounded by the superstructure”. A more gradual 
forming oxide layer on Cu (111) is presented in Figure 2 where one sees a set of images 
obtained at a constant potential for 10 minutes. The oxide free surface (-800 mV, t = 0), 
image 2a, undergoes a change involving lateral terrace growth and oxide island formation 
within the first two minutes (image 2b). The oxide layer during the next eight minutes 
(images 2c-2e) covered a larger portion of the surface where most of the terrace regions 
are covered by the UPD oxide. The pit-like structure is interpreted by Maurice et al. [15] 
as a “mixed layer containing metal and oxide species”. According to these authors the 
“pits” are not true topographic features. They argue that islands of adsorbate with lower 
apparent height are due to the lower conductivity of the oxide. Almost complete surface 
recovery was observed in the last image (2f) where a potential of — 1300 mV was applied 
to sample. 
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Summarizing the results of the ECSTM experiments it was concluded that an UPD 
oxide layer covers the Cu (111) surface in the potential range of (-S00 mV to -750 mV) 
in a pH 10 electrolyte. This layer is stable at potentials positive to -800 mV and is 
completely reduces at -1300 mV. Furthermore, the electrical conductivity of the Cu 
surface covered with the UPD oxygen adlayer is high enough as to allow high quality 
STM imaging. 

RDE results 

Figure 3 shows sets of anodic polarization curves obtained on Cu and Pt 
polycrystalline electrodes in 0.1 M Na 2 S 04 (pH = 5.5, Figure 3a) and in 0.1 M Na 2 S 04 
adjusted to pH 10 (Figure 3b) solutions. It is obvious that under these conditions the 
electrochemical behavior of Cu and Pt as substrates for the oxygen reduction reaction is 
identical and independent of pH. The diffusion-limited current densities measured on Cu 
and Pt are plotted as function of the square root of the disk rotation rate (to) in Figure 4. 
The slopes in Figure 4 are 78.1 x 10' 6 A s 1/2 rad' 1/2 cm' 2 for Cu and 77.7 x 10" 6 A s 1/2 rad' 1/2 
cm' 2 for Pt. The linear dependence of the current density on (0 I/2 is indicative of diffusion 
limited process that is described by the Levich equation, i L = 0.62 nFD 2n v~ V6 c^co U2 . 
Here k is the diffusion limited current density, n is the number of transferred electrons, F 
is Faraday’s constant, D the diffusion coefficient, v the kinematic viscosity, to the disc 
rotation rate and Coo the concentration of oxygen in the bulk solution. The experimentally 
obtained slopes on Cu and Pt are very close to the theoretically calculated slope (83.3 x 
10' 6 A s 1/2 rad' l/2 cm" 2 ) for a four electron (n = 4) reduction process for D = 1.93 x 10" 5 cm 2 
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s' 1 [16], v =1.01 x 10‘ 2 cmV [17], and c^-2.25 x 10" 7 mole cm" 3 [18] in naturally 
aerated solutions. Markovic et al. used a rotating ring-disk electrode configuration to 
study of oxygen reduction on Pt ( hkl) [19] in oxygen saturated 0.05 M H2SO4. They 
demonstrated that oxygen reduction on Pt occurred via the 4-electron mechanism and 
obtained a Levich slope of 463 x 10‘ 6 A s 1/2 rad" 1/2 cm" 2 which agrees well with the 
theoretically calculated value (466 x 10" 6 A s l/2 rad" 1/2 cm" 2 ) for n = 4. This slope is larger 
then the one obtained in the work presented herein by a factor of 5.9. This difference 

-7 

correlates with the difference in the concentration of oxygen in the electrolyte, 12.6 x 10" 
mole cm" 3 [17] for the oxygen - saturated solution and 2.25 x 10" 7 mole cm" 3 [18] for the 
naturally aerated solution. Similar comparisons exist with results of Vasquez et al. [13] 
who investigated oxygen reduction on a pure copper electrode in oxygen - saturated 
borate solution. In that work the authors determined n = 3.8 for the electron pathway of 
the reaction. This result was obtained under the assumption that D = 1.9 x 10"' cm s" , v 
=1 x 10" 2 cmV, and c x is 12.6 x 10" 7 mole cm" 3 . We have estimated the slope of the 
Levich plots presented in [13] to be 440 x 10" 6 A s 1/2 rad" 1/2 cm" 2 . This slope is close to 
that obtained by Markovic et al. [19] on Pt in 0.05 M H 2 S0 4 and differs from the one we 
obtained in the present study by a factor of 5.6. These results taken together lead us to 
conclude that a four-electron oxygen reduction mechanism takes place on copper in 
naturally aerated 0.1 M Na 2 S0 4 solution. Interestingly enough, we found no evidence of 
an influence of either oxygen concentration or the existence of the UPD oxygen adlayer 
on the reduction mechanism or the magnitude of the diffusion limited current density. 
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Figure 5 shows four sets of polarization curves obtained in 0.1 M Na 2 SC >4 solutions 
at different pH levels (3; 5.5; 7.5; 12). Each set contains polarization curves registered at 
prescribed disk rotation rate. The curves represented in each of the sets with solid lines 
were obtained while sweeping the potential from open circuit in cathodic direction. The 
dotted line curves were initiated at the onset of hydrogen evolution and the potential was 
swept in an anodic direction. When the experiment is performed in alkaline solution, 
there exists a wave within the potential range -600 mV to -400 mV. As seen in Figure 5, 
this wave does not exist in acidified or neutral solutions. According to the Pourbaix 
diagram [20], this wave can be attributed to the oxidation of Cu to Cu(I). The major 
difference in the anodic and cathodic sweeps is due to the initial chemical nature of the 
copper surface. In an anodic sweep the initial Cu surface is film-free regardless of the 
solution pH. In the cathodic sweep case at pH > 6 [20] the initial copper surface (i.e., at 
open circuit) is covered with a bulk oxide (CU 2 O). Comparison of the curves at different 
pH levels shows almost identical plateau-levels on the anodic and cathodic curves at pH 3 
and pH 5.5. The plateau-levels are different at pH 7.5 and pH 12 in that the cathodic 
curves systematically exceed the anodic ones. The difference in it (the plateau-levels) is 
summarized in Figure 6 where all the limiting currents are plotted as a function of pH and 
sweep direction. There is no difference in ii in the anodic and cathodic sweeps in 
acidified and neutral solutions (pH 1.5 to 5.5). In alkaline solution the current densities 
measured by a cathodic sweep exceed the ones obtained by anodic sweep by 50 to 150 
pA cm' 2 . This result is most likely due to a contribution of bulk oxide reduction to i L . 
According to [15], the reduction of Q 12 O in 0.1 M NaOH takes place between -600 mV 
and -1000 mV; i.e., within the range of our measurement of ii. Since bulk oxide 
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reduction is not under diffusion control its contribution is independent or only slightly 
dependent on the disk rotation rate. No influence of the UPD oxygen adlayer on the 
oxygen reduction reaction could be detected in any of our results. 

Figure 7 shows Levich plots for it as a function of rotation rate and pH. Figure 8 
summarizes the results for the slopes of the Levich plots shown in Figure 7. The average 
slopes of the lines obtained from the cathodic and anodic polarization curves are 83 ± 7 x 
10' 6 A-s l/2 rad' l/2 cm" 2 and 79 ± 3 x 10' 6 A-s 1/2 rad' l/2 cm" 2 respectively. This result is 
consistent with the theoretically predicted Levitch slope of 83.3 x 10 A-s 'rad cm" 
for a four-electron reaction (n = 4) in naturally aerated solution. The results are entirely 
consistent with those obtained by Markovic et al. on Pt [19] and therefore indicate the 
operation of a four-electron mechanism for oxygen reduction on Cu over the entire pH 
range examined. 

4 . Summary and Conclusions 

The reduction of oxygen on a polycrystalline copper electrode was investigated in 0.1 
M Na 2 S0 4 solutions at different pH levels. RDE and ECSTM techniques were used to 
examine the Cu surface morphology and the kinetics of the process. 

The ECSTM study showed a formation of a UPD oxygen adlayer on Cu(l 11) at 
potentials far more negative (-800 mV) then the potential of Cu(I) oxidation in a pH 10 
electrolyte. In agreement with the results of Maurice et al. [15] obtained in a pH 13 
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electrolyte, the evolution of that layer at a fixed potential was found to depend on the 
electrochemical potential and time. No influence of the adlayer on the oxygen reduction 
reaction kinetics was detected over the entire pH range examined. This indicates that the 
conductivity of this layer was high enough as to not impede the mass-transport limited 
rate of oxygen reduction. 

Results of RDE experiments show that a four - electron mechanism for oxygen 
reduction on Cu is operative over the entire pH range investigated. The initial stage of 
the surface (film - free Cu or bulk QhO) was found to slightly affect the diffusion- 
limited current densities measured at pH > 6. 

It is probably worth mentioning that we performed a similar set of experiments 
measuring i l in a thin - fluid layer cell under stagnant conditions and obtain the same 
value as those quoted above. These thin - layer cell experiments essentially maximize 
the electrolyte/vapor interfacial area with respect to the cell volume. All of the results 
taken together lead to the conclusion that the diffusion limited current density for oxygen 
reduction via a four - electron mechanism in a stagnant electrolyte is 20 - 30 |iAcm' . 
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Figure 1. STM micrograph showing morphological changes on Cu (111) in 0.1 M 
Na 2 S 04 adjusted to pH 10 solution as a function of increasing potential; (a) -1100 mV, 
(b) -1000 mV, (c) -900 mV, (d) -800 mV, (e) -750 mV, (f) -1300 mV 

Figure 2. STM micrograph showing morphological changes on Cu (111) at a constant 
potential of-800 mV in 0.1 M Na 2 S0 4 adjusted to pH 10 solution as a function of time, 
(a) 0 minute, (b) 2 minute, (c) 5 minute, (d) 8 minute, (e) lOminute and (f) -1300 mV. 

Figure 3. Anodic polarization curves obtained in: (A) and (B) 0.1 M Na 2 S 04 adjusted to 
pH 10 on Cu (solid lines) and Pt (dashed lines) polycrystalline disks at pre-described disk 
rotation rates. 

Figure 4. Levich plots for oxygen diffusion-limited current densities measured in 0.1 M 
Na 2 S0 4 and 0.1 M Na 2 S0 4 adjusted to pH 10 on Cu and Pt polycrystalline disks at pre¬ 
described disk rotation rates. 


Figure 5. Anodic (dashed lines) and cathodic (solid lines) polarization curves obtained in 
0.1 M Na 2 S0 4 as a function of pH (A, pH = 3; B, pH = 5.5; C, pH = 7.5; D, pH = 12) on 
copper polycrystalline disks at pre-described disk rotation rates. 
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Figure 6. Oxygen diffusion-limited current densities (anodic - dashed lines and cathodic 
- solid lines) measured in 0.1 M Na 2 S0 4 as a function of pH on copper polycrystalline 
disks at pre-described disk rotation rates. 


Figure 7. Levich plots for oxygen diffusion-limited current densities measured in 0.1 M 
Na 2 S0 4 as a function of pH on copper polycrystalline disks at pre-described disk rotation 
rates; A, i L measured during an anodic sweep; B, i L measured during a cathodic sweep. 


Figure 8. Summary of slopes of the Levich plots presented in Figure 7as a function of 
pH. 
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